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ABSTRACT 


A gray-level transformation method suitable for inter- 
active image enhancement is presented. It is shown that the 
well-known histogram equalization approach is a special case 
of this method. A technique for improving the uniformity of 
a histogram is also developed. Experimental results which 
illustrate the capabilities of both algorithms are described, 
Two proposals for implementing gray- level transformations 

in a real-time interactive image enhancement system are also 
presented. 
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CHAPTER I 


INTRODUCTION 

The O.VB& of research known as digital image processing 
can be broadly divided into two categories. The first cate- 
gory includes image processing for feature extraction and 
classification. The second category includes image restora- 
tion, enhancement and coding. 

The first category deals with techniques such as edge 
and curve detection and texture measurement. These are 
operations which do not necessarily output a new image but 
rather detect features within an image that can be analyzed, 
for example, by some typ of pattern recognition device. 

Image coding is primarily concerned with the techniques 
used in transmission and storage of images in a digital for- 
mat. In general, the objective of this type of work is to 
reduce the number of bits per picture element needed to accu- 
rately reproduce the image. In transmission, the objective 
is to reduce the bandwidth required in the communication 
channel. In the storage application of coding, the objec- 
tive is to reduce the amount of bulk storage required for an 

image . 

Image restoration is cor 'erned with modifying an image 
that has been degraded in some way, in an attempt to arrive 
at an estimate of the original image. This degradation 
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could take the form of noise injected over a communications 
channel, motion blur or lens aberrations in the imaging 
system. 

Image enhancement is an area of work concerned with 
manipulation of an image to extract details that might not 
be obvious in the original picture. Enhancement techniques 
can be broken down into two general classes: (1) transform 

domain methods, and (2) image domain methods. The first 
class includes techniques based on operations performed on 
a two-dimensional transform of an image such as, for example, 
the Fourier and Hadamard transforms. This type of processing 
is well-suited for removal of periodic noise, sharpening and 
smoothing operations. 

The second class of enhancement techniques are those 
that operate on the image itself. These techniques can be 
further broken dovm into position-variant and position-invari- 
ant methods. The position-variant methods are those that 
view the image as a two-dimensional array of points. Hence, 
the operation performed on a point depends on its position 
in the picture. Position-invariant methods of enhancement 
are the subject of this work. These are the methods that 
operate directly on a picture element, without regard to its 
position in the picture. One class of these methods is 
known as gray- level mapping. The objective of this approach 
is to find a function that will map the gray levels of an 
image into another set of levels to enhance some feature of 
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the image. As will be seen in the following chapters, this 
approach readily lends itself to implementation in an inter- 
active image enhancement system which requires no knowledge 
of enhancement techniques on the part of the operator . 

A brief review of the terminology and image formats 
used in this report will be useful in establishing the nota- 
tion used in subsequent chapters. A picture element will be 
referred to as a pixel. It is a single point in the picture 
and can assume any one of a number of quantized values. In 
the system used for this work, each pixel is represented by 
eight bits or 256 discrete values. The values can be con- 
sidered to span the range (0,255j or they can be normalized 
to a range of [0,1). At times it may also be advantageous 
to view the range of pixel values as {-0.5,0.5] . The repre- 
sents ion used depends on the particular application. 

The system used in the experiments to be described in 
subsequent chapters is shown in Figure 1. It consists of a 
Digital Equipment Corporation (D.E.C.) PDF 11/40 minicomputer 
with thirty- two thousand words of core memory. Standard 
peripherals interfaced to the system include two high density 
disk drives, an industry-compatible magnetic tape unit, an 
alphanumeric display terminal, a high speed paper tape 
reader/punch, a card reader, and a line printer. 

The image-inputting device consists of a standard 
closed-circuit television camera and a video sampling system. 
The sampling system consists of a hardware unit that monitors 
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the horizontal and vertical synchronization signals from the 
television camera to generate a series of pulses that ini- 
tialize an analog to digital converter to sample the video 

signal. It also consists of the software to drive the 

( 2 ) 

system. ' 

The video display system used for outputting the proc- 
essed results consists of a Model 6600 Display System manu- 
factured by Data Disc, Inc. This system generates a video 
signal with 480 by 640 individually addressable points. For 
the purposes of this work, the pictures were displayed in a 
480 by 512 format. The video display unit as it is now con- 
figured generates a 256 level monochrome signal, and a 256 
level color signal that drive standard black and white and 
color television monitors. 

The sampling unit and all of the software developed thus 
far are capable of handling three different picture sizes : 

128 X 128, 256 x. 256, and 512 x 512. These will be referred 
to as small, medium and large, or S, M, and L. 

Previous work done on gray level transformations will 
be discussed in the following chapters, as well as the algo- 
rithms that have been developed during the course of this 
investigation. The considerations relative to the inter- 
active nature of the processing will also be discussed. 


Superscript numbers in parentheses refer to those 
entries listed in the List of References. 



Finally, experimental results will be discussed and conclus 
ions will be drawn as to the value of this work and what 
might lie in the future for this type of image processing. 


CHAPTER II 


GRAY -LEVEL TRANSFORMATIONS 


Background 

The objective of the work in gray-level transformations 
is to find a function that will modify the gray levels of an 
image, in order to enhance some feature of the image. Let- 
ting X represent the gray levels of an image to be enhanced, 
the transformations of interest are of the form 

y = T(x) (1) 

which map the levels x into corresponding levels y to obtain 
an image which has been enhanced in some manner . It is 
assumed in the following discussion that the levels have 
been normalized to the interval [0,1], where zero represents 
black and one represents white in the gray scale. 

Two conditions must be imposed on the transformation 
function T(x) to obtain meaningful resul:.>’. First, T(x) 
must be monotonically increasing in the sense that 

T(x^) < TCx^) < X 2 >. (2) 

The implication of this constraint is that if pixel Xj^ is 
lighter than pixel X 2 in the original image, then pixel x^ 
will also be lighter than pixel x^ in the enhanced image. 
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If this condition is not met, the resulting picture could be 
meaningless . 

The second condition imposed on the transformation is 

that 

0 < T(x) <1 (0 < X < 1). (3) 

This serves to maintain a mapping that will be consistent 
with the allowed range of gray levels. 

The variable x is a random variable. For each image to 
be enhanced, x has a probability density function or p.d.f. 
denoted p (x) . In the enhanced image, the levels will be 
denoted by y with p.d.f. Py(y) . For a discrete or digitized 
image, the p.d.f. is a discrete function known as a histogram 
which specifies the number of occurrences of each gray level 
within an image. For simplicity, the following development 
is presented in continuous mathematics. The discrete equiv- 
alents can be obtained without difficulty. 

Given a transformation function T(x) and a density p (x) , 
statistical theory provides a means for estimating the p.d.f. 
of the random variable, y, given that y = T(x). The p.d.f. 
of the variable y in Equation (1), page 7, can be obtained 
using the relation 


- ... S _ 

Py^^i^ |T’ (Xj^) I I T' I 


(4) 
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where y. = T(x.); j = 1, n; and T'(x) denotes the derivative 
^ J 

( 3 ) 

of T with respect to x. Since one of the constraints on 
the transformation T is that it be strictly increasing, there 
will be only one solution to Equation (1) , page 7 . In addi- 
tion, the absolute value sign may be dropped since T'(x) will 
always be positive because of this condition. Hence, Equation 
(4) , page 8, becomes 


Py(y) = 


P^(x) 
r (x) 


P^(X) 




= p- 




(5) 


where x = T ^(y). Thus, the p.d.f, of the enhanced image 
can be estimated knowing only the original p.d.f. and the 
tr an s forma t i on f unc t i on . 

It is noted that knowledge of the function p (x) and 
Py(y) yields information about the physical characteristics 
of the images. For instance, the average brightness, or 
mean, of the original image is given by 


X =f X p (x)dx, 
JQ ^ 


( 6 ) 


while the variance, 


‘’x ' Ig (X-J)2p^(x)dx, 


(7) 


is a measure of the amount of contrast in the image. The 
brightness and contrast parameters are of fundamental impor- 
tance in describing the characteristics of an image. They 
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also provide a basis for comparing images mathematically 
using parameters that can be interpreted visually. 

Histogram Equalization 

A very popular approach to the enhancement problem is 
to generate a uniform p.d.f. The objective in this case is 
to redistribute the gray levels in a picture such that each 
gray level within the p.d.f. contains the same number of 
pixels. The intention is to optimize the use of the avail- 
able levels. This technique is known as histogram equaliza- 
tion or distribution linerization . 

The transformation used is of the form 


y = T(x) = [ p (s)ds (0 < X < 1) , (8) 

io 

where s is a dummy variable of integration. This function 
is known as the cumulative distribution function or c.d.f. 
of X. From Equation (8) , it follows that 


dx _ 1 

dy Pjj(x) ’ 

therefore, Equation (5), page 9, reduces to 

Py(y) =1 (0 < y < 1) 


(9) 


( 10 ) 


which is the desired uniform density. 

A heuristic argument can be presented in favor of histo- 
gram equalization. It is not unreasonable to assume that 
the quality of an image will be generally improved if its 
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pixels are evenly distributed over the entire available gray 
scale. All of the significant details present in the origi- 
nal image are generally preserved, while the subtle details 
not obvious in the original image are enhanced. This can 
best bo illustrated by example. 

Figure 2 is the image of a woman's face with a shadow 
cast over one side. Figure 3(a) illustrates the histogram 
of the image before processing. The contrast of the image 
is rather poor, as would be expected from the narrow 
gray-level range spanned by its histogram. The c.d.f. of 
this histogram appears in Figure 4. The reference line is 
to compare the actual c.d.f. with the c.d.f, of a uniform 
density which is the curve y = x. 

The image of Figure 2 was transformed using the curve 
in Figure 4, The result appears in Figure 5, with its corre- 
sponding histogram in Figure 3(b). The transformation has 
tended to flatten and expand the histogram. The gaps in the 
resulting histogram are due to the fact that the pixels 
assume only discrete values. The c.d.f. of the resulting 
histogram appears in Figure 6. Note that it more closely 
approximates the reference line. The histogram indicates 
that the contrast of the image has been increased and Figure 
5 confirms this point. Note that there is considerably more 
detail visible in the enhanced image than in the original. 

For most images, this has proven to be the case. 
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Figure 2. Original image of a woman's face with sidelighting. 

Figure 3. Histogram of Figure 2 (a) before processing and 
(b) after equalization. 

Figure 4. Cumulative distribution function for Figure 2. 
Figure 5. Histogram-equalized version of Figure 2. 

Figure 6. Cumulative distribution function for Figure 5. 
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Contrast Enhancement Methods 

Another type of gray level mapping is referred to as 

(5 7) 

contrast enhancement or contrast stretching. * The 
objective of this type of processing is to maiiipulate the 
spread of one or more areas of the density function. Some 
example transformation curves appear in Figure 7. Figure 7 
(a) shows a transformation function that could be used on an 
image where all or most of the pixels have values on the 
interval ta,b). This transformation then spreads or 
stretches the density to cover the entire [0,1] range. 

Figure 7(b) shows a curve that could be used to spread the 
dark region of a picture, while Figure 7(c) illustrates a 
curve that will operate similarly on the light areas of a 
picture. Finally, Figure 7(d) shows a curve that would tend 
to spread out the middle range of a picture. 

One algorithm has been proposed that generates a func- 

/ o\ 

tion similar to that in Figure 7(d). The function used 

is of the form 

y = T(x) = f“^(tan"^(c*tan(f (x)))) (x 0 or 1) , (11) 

and 

T(0) = 0, T(l) - 1, (12) 

where f is a function that maps the pixel spectrum [0,1] onto 
the domain of the tangent function [ -tt /2 ,ti / 2 ] . The constant 
c determines the degree of enhancement. 


Figure 7. Transforrnation functions for (a) increasing the 

dynamic range of a picture, (b) spreading the daric area 
(c) spreading the light area, and (d) spreading the 
middle range. 
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The tangent function maps its domain onto the infinite 
real line. The positive constant c multiplies the points on 
the real line to either spread them out or compress them. 

Then the inverse tangent and inverse f are taken to map the 
points back to the interval [0,1]. Hence, the value of the 
constant c directly controls the variance or spread of the 
density function of the enhanced image. 

The function f, as it was propose!, maps the mean of 
the density to the raid-point of the interval [-ti/2,ti/2] which 
is zero. Similarly, the inverse f maps the mid-point back 
to the mean of the original image as explained in Reference 
8. The inverse of f can be replaced by another function 
similar in nature to the inverse of f. This function can be 
used to map the mid-point of the tangent domain to a mean 
specified by an operator, rather than the original mean. 

This gives the user control over the mean as well as the 
variance of the result. For the following examples, it is 
assumed that the old mean and the new mean are both 0.5. 

Figure 8(a) is the transformation function defined in 
Equation (11), page 14, that was generated with c = 1.0, 
which yields the original image. Figure 8(b) was generated 
with c = 0.2. Assuming that most of the pixel values lie in 
the mid-range of the specui.am, this transformation will tend 
to decrease the contrast of the picture. Figure 8(c) shows the 
function that was generated with c = 5.0. This transformation 
will tend to increase the contrast of the picture. Note 




'8 

that each of these curves has odd symmetry about the point 
(0.5, 0.5). This is due to the fact that the tangent func- 
tion has odd symmetry about the origin. 

Figure 9 is the result of contrast enhancement as 
proposed in Reference 8. The original of this image appears 
in Figure 2, page 12 . Fairly good results have been obtained 
with this algorithm in most cases. This approach gives the 
user direct control over the contrast of a picture through 
the parameter c, and the brightness through the function f. 

Drawbacks of Previous Methods 

Two methods of gray-level transformations have been pre- 
sented with examples to show their effectiveness. Arguments 
are given below to show the inadequacies of each method in 
terms of implementation in an interactive environment and in 
terms of the algorithm itself. 

The major drawback of histogram equalization is that 
the algorithm is fixed and, therefore, the user has no con- 
trol over the process. It will be shown below, however, 
that histogram equalization can be used as an intermediate 
step in an algorithm that provides results superior to 
either of the methods described previously. Further, it will 
be shown that histogram equalization is a special case of 
this method. 

The contrast enhancement approach proposed in Reference 
8 lends itself well to implementation rn an interactive 
environment, requiring only two parameters to control it. 
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Figure 9. Contrast enhanced version of Figure 2 
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The algorithm is fairly straight-forward and does not involve 
much computation. However, it has one major flaw. As 
pointed out above, the function as defined in Equation (11), 
page 14, is symmetric about the mean of the original image 
as was illustrated in Figures 8(b) and 8(c), page 17. This 
implies that the density of the image to be enhanced should 
be symmetric about its mean, a condition which is seldom 
satisfied in practice. This algorithm is certainly suitable 
for manipulating the brightness and contrast of pictures, 
however, this is the extent of its usefulness. It lacks the 
flexibility required to manipulate an unsymmetric density ir 
an unsymmetric manner. 


Direct Histogram S p ecification 

The gray-level mapping method developed below, like 
those presented above, is based on transforming the prob- 
ability density function of an image to be enhanced. The 
process is different, however, in that a transformation is 
developed that will produce an image with a density function 
py(y) which has been interactively specified by an operator. 

Using Equation (8), page 10, the gray levels x of an 
image can be transformed to a new set of levels z, that is, 


z = H(x) = [ p (s)ds (0 < X < 1) . (13) 

From the above discussion it is known that p„(z) will be a 

/it 

uniform density function. If the inverse of Equation (13), 
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H"^(z) is applied to the new density, the original image 
with density p (x) is obtained. It is noted, however, that 
if Equation (8), page 10, is applied to a different density 
Py(y) > the result will be 

2 = G(y) = p (s)ds (0 < y < 1). (14) 

jq y 

Although G(y) is in general different from H(x), is 

the same if either Equation (13), page 20, or Equation (14) 
is used. If the inverse function G”^(z) is applied to the z 
levels, the result will be a set of y levels with the speci- 
fied density p^(y) . In other words, if the z levels have a 

uniform density, the desired p (y) can be obtained using the 

-1 -1 
inverse mapping G (z) . It is noted that both H(x) and G (z) 

must satisfy the condition of monotonicity as set forth in 
Equation (2), page 7. Therefore, if th:j original image is 
first histogram-equalized, and the new (unifonn) levels are 
inverse-mapped using the function G ^(z) , the result will be 
an image whose gray levels have the desired density py(y) • 

The procedure is illustrated in Figure 10. Figure 10(a) 
shows the original density function Pjj.(>^) which is computed 
from the image to be enhanced. Figure 10(b) shows the cumu- 
lative distribution function H(x) which is used to map the 
variable x into the variable z. Figure 10(c) illustrates 
the result of this transformation. Figure 10(d) is the c.d.f. 
of the desired density, G(y) , the inverse of which is used 
to map the variable z back into the desired density p^(y). 
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Figure 10(e) illustrates the result of this transformation 
which is the desired density. 

This procedure can be expressed in the form of an equa- 
tion that specifies the transformation from the given density 

p (x) to the desired density p (y) . The relation is given by 
X y 

y=G'^[H(x)l (15) 

where H(x) is the c.d.f. of x and G(y) is the c.d.f. of y. 

It is noted that histogram equalization is a special case of 

“3 (9) 

this technique obtained by letting G [H(x)] = H(x) . 

This method of transformation from one density to 
another is well known to statisticians. The problem in using 
this technique with continuous functions is that for most 
functions, the inverse of the c.d.f. is difficult if not 
impossible to evaluate analytically. For discrete functions, 
however, the inverse can be approximated without difficulty. 
Since the emphasis of this investigation is on discrete 
images, it follows that Equation (15) can be used as a tool 
for image enhancement, provided that a desired hittogram has 
been specified. As will be shown in the following chapter, 
the problem of specifying a histogram can be formulated in 
an interactive manner, thus giving the operator control over 
the enhancement process. 

As an example of the superiority of this method over 
histogram equalization, consider Figure 11 which shows the 
image of a quarter taken inside a poorly illuminated box. 
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Fir.uro 11. Original image of a quarter in 
illuminated box. 
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Figure 12 is the result of histogram equalization. Figure 

13 illustrates the result of contrast enhancement as pro- 
posed in Reference 8. Figure 14 is the result of processing 
with a method that uses the transformation given in Equation 
(15) , page 23 . The method used to specify the resultant 
density will be described in Chapter III. Note that Figure 

14 provides much more information about the object in the 
image than do Figures 12 or 13. The histograms for Figures 
11, 12, 13, and 14 appear in Figures 15(a), 15(b), 16(a), 
and 16(b), respectively. It is noted that the histograms of 
the three enhanced pictures do not differ greatly, however, 
the amount of visual information in the picture does . 

Improving the Uniform Density 

The success of the above algorithm in approximating a 
specified density relies on the fact that a uniform histogram 
is achieved in the process of histogram equalization. Due 
to the fact that the pixels assume only discrete values, 
however, the result of the equalization process is generally 
only an approximation to a uniform histogram. Figures 3(a) 
and 3(b), page 12, illustrate a good example of the non-uni- 
formity that is achieved in practice. If the original density 
has even fewer gray levels present, the approximation gets 
worse as is illustrated in Figures 15(a) and 15(b). The 
root mean square (R.M.S.) error between each of these equal- 
ized histograms and a true uniform histogram is ,0063 and 
.031, respectively. 
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Figure 12. Histogram-equalized version of Figure 11. 

Figure 13. Contrast enhanced version of Figure 11. 

Figure 14. Image of Figure 11 after enhancement by histogram 
/ specification. 

Figure 15. Histograms of (a) Figure 11 and (b) Figure 12. 

Figure 16. Histograms of (a) Figure 13 and (b) Figure 14. 



Figure 12 
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Figure 15 


Figure 16 
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Due to the discrete nature of the pixel values, even if 
a uniform histogram is obtained, it will not necessarily 
yield exactly the histogram that was specified. Theoreti- 
cally, however, this is the best that can be achieved. 

Hence, a uniform histogram can be used as a reference in 
determining how close the specified and resultant histograms 
actually are. For example, Figure 17 illustrates a histogram 
that v;as specified using the method to be outlined in Chapter 
III. Figure 18 shows the result of mapping a perfect, uni- 
form histogram to the histogram shown in Figure 17. The 
R.M.S. error between the specified and resultant histograms 
was computed to be .00'' 7, Figure 19 illustrates the result 
of mapping the histogram of the image in Figure 2, page 12 , 
to the histogram specified in Figure 17. The R.M.S. error 
was computed to be .006 which is a 243.8 percent increase 
over the error obtained vi;ith the uniform histogram. 

An algorithm has been developed that can be used to 
preprocess an image in order to improve the uniformity of 
Che histogram resulting from the equalization procedure. 

The algorithm was adapted from some resu] ts in texture- 
analysis . This algorithm differs from those previously 

described in that it considers a pixel and a small neighbor- 
hood around it, rather than just the pixel itself, when 
determining the new value of the pixel. 

The image is assumed to consist of a two-dimensional 
array A with elements A(i,j). A new image B is computed 



le of a hisrop,raiii specified by a usei’. 


nj; of a uniform hisLoj-ram to that in 



m’ of thc' image in Fig.uro 2 to the 
Figure 17. 
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from A using the following relation: 

1 1 

= I I A(i-m.j-n) (16) 

m=-l n=~l 

assuming that the pixel values are on the interval [0,255). 
This yields an image with a resolution of one part in 2,304 
rather than one part in 256. This is effectively ai^ increase 
from eight bits to approximately 11.2 bits. This new image 
B with a 2,304 point density function can then be histogram- 
equalized and quantized hack into 256 levels. 

The significance of this algorithm is that it can dif- 
ferentiate between two different pixels with the same value 
that appear in the image in different contexts. For example, 
suppose that pixels and X 2 have the same value, but pixel 
x^ is situated in a very dark region of the picture while 
pixel X 2 is in a lighter region. This algoritlam, rather 
than map both x^ and X 2 to the same value, will map to a 
dar'ker value than X 2 • 

The reason that histogram equalization does not yield a 
uniform density is that for any gray level, all of the pixels 
with that value will be mapped to the same new value. This 
algorithm, however, provides a way of mapping contextually 
different points with the same value to different values. 

This allows large groups of points with the same value to be 
broken dovm into smaller groups with different values in a 
meaningful manner. 
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One of the advantages of this preprocessing technique 
is that it tends to average out the noise in an image. 

Another advantage is tnat successive iterations through the 
algorithm tend to produce better approximations to the 
de.sired uniform histogram. Hence, the histogram obtained 
in the second step of the algorithm will more closely approx- 
imate the histogram that was specified. Figure 20 illus- 
trates the histogram of Figure 2, page 12, after one itera- 
tion through the preprocessing algorithm. The R.Ii.S. error 
of this histogram V7ith respect to a uniform histogram is 
.0008A. The use of this preprocessed image when input to 
the histogram specification algorithm results in the histo- 
gram shown in Figure 21. The R.M.S. error for this histo- 
gram was computed to be .00195, which is only a 12.2 percent 
increase over the reference. Figures 22 and 23 illustrate a 
second iteration through the preprocessor. The R.M.S, error 
of Figure 22 with respect to the uniform histogram in this 
case has decreased to .00029. The error in the resultant 
histogram has decreased to .00186, which is only 7,0 percent 
above the reference. Note that both iterations have reduced 
the R.M.S. error at ’'Oth the first and second steps of the 
specification algorithm. Experience has shown that after 
two iterations the histogram docs not improve significantly. 

There is a problem with the preprocessing method des- 
cribed above. Like neighborViood averaging, it tends to blur 
the picture. Figure 24, for example, illustrates a picture 
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Figure 20. Histogram of the image in Figure 2 after one 
iteration through the preprocessor. 

Figure 21. Mapping of the histogram in Figure 20 to the 
histogram in Figure 17. 

Figure 22. Histogram of the image in Figure 2 after two 
iterations through the preprocessor. 

Figure 23. Mapping of the histogram in Figure 22 to the 
histogram in Figure 17. 



Figure 20 



Figure 22 





Figure 2k. Original image of a woman's profile. 
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before processing of any type. Figure 25 shows the same 
picture after histogram equalization. Figures 26 and 27 
show the results of one and two iterations of preprocessing, 
respectively. Note that the detail has been blurred to 
some extent. 

A rather heuristic method has been developed for elimi- 
nating the blur without significantly affecting the general 
shape of the histogram. Let C(i,j) represent nine tir.3s 
the value of A(i,j). (It takes nine points from A to yield 
one point in B.) After the point B(iJ) is computed as in 
Equation (16), page 30, the difference between this point 
and C(i,j) is computed, A histogram of the absolute value 
of these diffe'^ences is then computed. The points where the 
difference is zero are disregarded. Since the reason that 
the picture blurs is that the computed value of B(i,j) dif- 
fers too much from C(i,j), it was decided to use the mean 
of this histogram as a threshold on either side of the 
point C(i.j). If B(i,j) lies v/ithin these limits, the point 
is not changed. If it lies outside of these limits, then 
the point is reassigned a value inside the limits. An 
attempt is made to distribute these values uniformly through- 
out the interval so as to aid in achieving a uniform density. 

Figures 28 and 29 illustrate the result of processing 
Figure 24, page 34, through one and two iterations of the 
preprocessor with thresholding. Note that there is no loss 
of detail with respect to the histogram-equalized image of 
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Figure 25. Histogram-equalized version of Figure 24. 

Figure 26. Image of Figure 24 after one iteration of the 
preprocessor . 

Figure 27. Image of Figure 24 after two iterations of the 
preprocessor . 

Figure 28. Image of Figure 24 after one iteration of the 
preprocessor with thresholding. 

Figure 29. Image of Figure 24 after tv;o iterations of the 
preprocessor with thresholding. 




Figure 25 


Figure 26 


Figure 27 


Figure 28 


Figure 29 
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Figure 25, page 36- Figures 30(a), 30(b), 31(a), and 31(b) 
illustrate the histograms of Figures 26, 27, 28, and 29, 
page 36. The R.M.S. errors were calculated to be .001, 
.00028, .00095, and .00048. Note that processing with and 
without thresholding makes very little difference in the 
histograms that were generated. It is obvious, however, 
that there is a very distinct difference in the quality of 
the images. 

In summary, an algorithm has been developed for obtain- 
ing a specified histogram from an arbitrary histogram using 
histogram equalization as an intermediate step. In general, 
the accuracy of the algorithm in achieving the desired 
histogram depends on the success of the histogram equaliza- 
tion in obtaining the desired uniformity. Hence, a preproc- 
essing technique has been implemented that acts to modify 
the image such that a better approximation to a uniform 
histogram can be achieved. Further results using this 
enhancement approach are presented in Chapter IV. 


Figure 30. Histograms of (a) Figure 26 and (b) Figure 27 


(b) 



(a) " 

Figure 31. Histograms of (a) Figure 28 and (b) Figure 29. 



CHAPTER III 


INTERACTIVE SPECIFICATION OF HISTOGRAMS 

As was stated in Chapters I and II, the principal objec- 
tive of this investigation was to develop an algorithm that 
could be implemented interactively. The reasoning behind 
this objective is that the quality of a picture is difficult 
to quantify. The computer can manipulate an image, but only 
a human with an intuitive knowledge of what is required can 
make the decision as to whether or not the computer has pro- 
duced a useful result. 

With the enhancement algorithm developed in Chapter II, 
all that is required is the specification of a desired histo- 
gram. The algorithm then computes a transformation function 
based on the histogram specified by the user and the histo- 
gram of the original picture. There are two basic categories 
of users that must be considered when devising a method for 
specifying histograms. The first category includes users 
who are familiar with statistical theory and image processing 
techniques and are interested in the theoretical aspects of 
the algorithm. The second category includes users who are 
only interested in the final results of the enhancement 
system. Both categories of users are considered in the 
following pages. 
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Very sophisticated schemes for specifying histograms 
can be developed when ic is known that the user has some 
degree of expertise in the field. One of the simplest 
methods would be to specify a histogram that approximates a 
normal density which requires only two parameters to define 
it. More complex schemes could include piece-wise linear 
histograms V7ith many segments. A procedure can be imple- 
mented v;hich will allow the user to experiment and determine 
what types of histograms produce certain features in an 
image . 

One of the more interesting methods of specifying a 
histogram involves merely sketching the desired curve. This 
curve can then be input to the computer in several different 
ways. A light pen could be used in conjunction with a 
graphics terminal for such an application. If a digitizing 
television camera is available, the curve can be sketched on 
paper and then digitized using some type of edge-following 
algorithm. Another method would be to use a joystick or 
trackball to move a cursor along the curve and digitize the 
X-Y position of the cursor. An example using this type of 
histogram input will be presented in Chapter IV. A fourth 
alternative would be to use a graphics tablet to sketch the 
curve . 

The possibilities for this type of specification are 
unlimited. The advantage is that the user can specify any 
shape of histogram that is required, without being limited 


by the structure of a histogram specification algorithm. 

From this standpoint, the user can learn what general shape 
of histogram yields a particular characteristic in a picture. 
An algorithm that simplifies the specification problem can 
then be formulated from this information. 

The methods mentioned above offer a wide range of 
options in the histogram specification process. The dravj- 
back is that they require knowledge of statistics and image 
processing techniques in order to be effectively u.sed. The 
untrained user may find it very difficult to extract usable 
output from such a system. A method of specifying a curve 
that is relatively simple, but effective, must be developed 
for this class of user. 

The first step in devising an algorithm that generates 
histograms is to detc nine the number of parameters that are 
required to adequately specify the major characteristics of 
a histogram. Two parameters that are obviously required are 
the mean and variance. As stated previously, the mean deter- 
mines the average brightness level of the picture, while the 
variance is a measure of the contrast of the picture, In 
terms of the histogram itself, the variance is a measure of 
the spread of the histogram about the mean. 

In the terminology of statistics, the mean and variance 
are referred to as the fir.st moment about the origin and the 
second moment about the moan, re.spectively . They are the 
only parameters required to completely specify the familiar 
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normal density. Hov/ever, the normal density is symmetric 
about the mean and as stated in Chapter II, a symmetric 
histogram is a rare case in practice. Something is needed 
to quantify the asymmetry of the histogram. 

The third moment about the mean of the histogram can be 
effectively used to express the asymmetry of a histogram. 

The third moment is given by the following relation: 


m 


3 


fl 
• 0 


(x-x)^ 


p^(x)dx. 


(17) 


There are several v;ays to measure the asymmetry or slcewness 
of the histogram. One is called the moment coefficient of 
ske\'7ness and is given by the relation 


m. 


^3 


(cjp3/2 


(IB) 


where a is the variance as defined in Equation (7) , page 9 . 
This coefficient v;ill be used throughout the following dis- 
cussion. There are two other coefficients of skev 7 ness 
kno^^n as Pearson's first and second coefficient of skewness. 
(11. 12) defined as follows: 


3 ( x-Mo) 
1 ■ 


(19) 


where Mo is the mode of the histogram and 
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Sk^ = (20) 

X 

where Md is the median of the histogram. 

For perfectly symmetrical histograms, all three of 
these coefficients are zero. If the coefficient is negative, 
the histogram is said to have negative skevmess or to be 
skev7ed to the left. This occurs when tlic liistogram ha.s a 
longer tail to the left of the central maximum than to the 
right, as in Figure 32(a). Figure 32(b) shows a histogram 
that is skewed to the right. In this case the skev7ness 
coefficient is positive, 

Most image histograms found in practice can bo ade- 
quately described statistically with the addition of the 
fourth moment about the mean (also called kurtosis), which 
is given by the relation 

m^ = (x-x)^ p^(x)dx. (21) 


As will be seen below, the kurtosis of a histogram can be 
used to describe its peakedness about the mean. The coef- 
ficient of kurtosis that quantifies this peakedness is given 
by the relation 



( 22 ) 


For the normal density, this coefficient is three. For this 
reason, the coefficient of kurtosis is sometimes defined as 




Using this relation, histograms can be separated into three 
categories. If the coefficient is negative, the histogram 
is termed platykurtic. An example of this appears in Figure 
33(a). If the coefficient is zero, the histogram is meso- 
kurtic and approximates the normal density. This corresponds 
to Figure 33(b). Figure 33(c) represents a leptokurtic histo- 
gram V7hich has a positive coefficient of kurtosis. 

The above development yields four parameters that can 
be used to describe a histogram. They give a method for 
describing the brightness and coiitrast of the image which are 
qualities that are visible. They also describe the asymmetry 
and peakedne.ss of the histogram which, although not readily 
interpretable in terms of visual quality, relate to the 
physical description of the histogram. Using the four sta- 
tistical parameters described above as a starting point, a 
method for specifying histograms interactively can be formu- 
lated. The basic approach is based on controlling the speci- 
fied histogram by manipulations which are related to these 
parameters . 

Foui parameters are used in the histogram specification 
method described below. The first parameter used is denoted 
by m. This specifies a point in the pixel spectrum about 
which the histogram will be generated. The second parameter 
is denoted by h and specifies the height of the histogram at 






the point m. The third and fourth parameters are referred 
to as the right spread (6^^) and the loft spread (0^) about 
m. It is noted that m is basically related to the mean of 
the histogram and h to its kurtosis. The other two para- 
meters control spread (variance) and symmetry. Although the 
number of parameters could easily be increased, experiments 
have indicated that more than four parameters are quite dif- 
ficult to control interactively, particularly when enhance- 
ment speed is an important consideration. 

There are two basic cases that illustrate the manner in 
which the histograms are specified. Figures 34(a) and 34(b) 
illustrate the approach for decreasing and increasing the 
spread about m. Both of these cases involve specifying the 
point m which can be located anyv?here on the interval [0,1], 
Then the height is chosen at that point and may assume any 
value greater than or equal to zero. 

For the case of decreasing spread, the parameters and 
specify the angle from vertical, as shov?n in Figure 34(a). 
These parameters can assume values from zero degrees to 90 
degrees. Point j moves along the line segment that connects 
(0,1) and (m,0) as 0j^ varies, while point k moves between 
(1,1) and (m,0) as 0^^ varies. For the case of increasing 
spread, 0j^ and 0^^ specify the points along the vertical axis 
at which the points j and k are fixed. The lower bound for 
these points is zero and theoretically there is no upper 
bound . 
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The four parameters are Sj ecified such that the left 
and right sides of the histogram can be manipulated inde- 
pendently. Once the parameters have been specified, a 
piece-wise linear curve can be constructed. The curve is 
then normalized to unity area v/hile preserving the shape of 
the histogram. Figure 35 illustrates four curves that were 
generated with this specification method. Figures 35(a) and 
35(b) illustrate positive and negative skewness as in Figures 
32(a) and 32(b), page 45. Figures 35(c) and 35(d) illustrate 
platykurtic and leptokurtic histograms as in Figures 33(a) 
and 33 (b) , page 47 . 

The above approach was implemented on a D.E.C. PDF 11/40 
minicomputer, using four system potentiometers which are 
interfaced to the processor. Each potentiometer controls 
one of the four parameters described above by varying a volt- 
age input to cin analog to digital converter. This voltage 
is then sampled and expressed as a number betv^een -i-1 and -1. 

If this algorithm were to be implemented in a real time image 
enhancement system, the potentiometers could be replaced by 
a joystick with four degrees of freedom as illustrated in 
Figure 36. This would be more convenient than the potentio- 
meters, and an operator with little or no knowledge of enhance- 
ment methods or the structure of the algorithm can rapidly 
learn to manipulate the joystick in order to achieve fast 
enhancement results. 
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(a) 


(b) 





(c) 


(d) 


Figure 35. Example histograms generated with the method 
described in Chapter III. 
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CHAPTER IV 


EXPERIMENTAL RESULTS 

This chapter presents several results of processing 
with the direct histogram specification algorithm developed 
in Chapters II and III. The pictures presented in the fol- 
lowing pages, as well as those of previous chapters, are 
all 256 X 256 or medium-sixed pictures. They are displayed 
in a 480 x 512 format, v.’hich constitutes a standard frame 
in a 525-line television monitor. The pictures are dis- 
played with 256 gray levels with no scaling or translation 
of the data before display. 

For each experiment, three pictures will be presented: 
(1) the original image, (2) the histogram-equalised image, 
and (3) the image that results from the histogram specifica- 
tion algorithm. Each of the pictures resulting from the 
specification algorithm was generated using two iterations 
of the nreprocessing technique described in Chapter II. The 
histograms of the original image, the histogram-equalized 
image, the desired image and the resulting image are also 
presented. They will be referred to as parts (d) , (e) , (f) 

and (g) , respectively. 

The four parameters used to generate the desired histo- 
gram of each picture (see Chapter III) arc listed in Table I. 
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TABLE I 

TABLE OF HISTOGRAM PARAMETERS USED TO 
GENERATE THE RESULTS OF CHAPTER IV 


Figure 

Number 

m 

0 i iTi < 1 

0 ^ 

G 

h 

0 i h 

37 

.350 

10.67" 

14.00 

0.32 

38 

.245 

4.91* 

7.44“ 

7.02 

39 

.630 

2.42* 

6.98* 

4.19 

40 

.770 

80.10 

•k 

9.21 

1.07 

43 

.705 

3.89* 

78.30 

0.78 

44 

.100 

49.50 

38.70 

0.40 


^Values of G, and 0^ denoted by an asterisk are to be 
interpreted as displacements along the vertical axis. (See 
Figure 34, page 49.) Otihcr values of 9^ and 0j^ are angles 
from vertical in degrees. 
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The values of 0, and Gr, that are denoted with an asterisk are 

la K 

to be interpreted as dispL- oments along the vertical axis. 
(See Figure 34, page A9.) The other values are to be inter- 
preted as angles from vertical in degrees. 

Result No. 1 

Figure 37(a) is the same quarter shown in Figure 11, page 
24- It is repeated here to illustrate the use of the preproc- 
essor. Figure 37(b) repeats the histogram equalization of 
the original picture. Figure 37(c) is the result of process- 
ing the preprocessed picture with the direct histogram speci- 
fication method. Note that the v;ord "LIBERTY" is discernible, 
whereas in the equalization case, it is almost completely 
obliterated. Figures 37(d) through 37(g) illustrate the 
histograms as outlined above. 

Result No. ?. 

Figure 38(a) is the image of a dark room as viewed from 
the doorway. The striped shirt of a subject is just barely 
visible. Histogram equalization produces the image in Figure 
38(b). Note that the contrast has been increased. Note also 
that the subject at the table has been enhanced, v/ith the 
chairs at the table partly visible. Further processing with 
the specification algorithm yields the image in Figure 38(c). 
Note that the chairs have been enhanced further and the 
details of the objects on the table have also been brought 
out. The histograms for this example are illustrated in 
Figures 38(d) through 38(g). 
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Result No. 3 

Figure 39 illustrates another view of the same room 
illustrated in Figure 38. In the original, Figure 39(a), 
there is not much contrast. Histogram equalization has 
increased the contrast, as illustrated in Figure 39(b), but 
has failed to bring out the detail in the room. Figure 39(c) 
illustrates the result of processing with the histogram 
specification algorithm. Note that the chairs in the room 
and some objects on the tables have been brought out in this 
example. Figures 39(d) through 39(g) illustrate the histo- 
grams pertinent to this example. 

Result No _4 

Figure A0(a) is a picture of a truck taken against a 
very light background. The brightness of the background has 
forced the dark details of the truck to have relatively very 
low contrast. Histogram equalization, in this case, served 
only to increase the contrast between the truck and the back- 
ground, as illustrated in Figure 40(b). Processing with the 
histogram specification algorithm produced the image in 
Figure 40(c) . The reason for the improvement in this picture 
is that the specification algorithm allowed the user to con- 
fine the background pixels to a very small portion of the 
spectrum, while spreading the dark pixels of the truck over 
the majority of the spectrum. Histogram equalization, on 
the othei* hand, spread the background pixels over most of 
the spectrum, leaving the important information in a very 
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Figure 40. Results of processing the image of a truck taken 
against a light background. 
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small region. This is illustrated in the histograms in 
Figures 40(d) through 40(g). 

It was mentioned in Chapter III that one way of speci- 
fying a histogram is to digitize a curve that has been 
sketched. The curve in Figure 41 was digitized using a 
television camera and a joystick. The program for doing 
this is listed in the appendix. Figure 42(a) illustrates 
the result of specifying the histogram in Figure 41. Figures 
42(b) through '(2(e) illustrate the relevant histograms in 
the same order as desci'ibcd above. 

Result No. 5 

The above examples illustrate images vjhose detail was 
enhanced by the use of the histogram specification algorithm. 
There are some cases, however, vjhere histogram equalization 
produces an acceptable image. Figure 43 illustrates mre of 
these examples. Figure 43(a) illustrates the image of the 
same truck as in Figure 40(a). This image, however, was 
taken against a dark background with very poor lighting. 
Histogram equalization succeeded in bringing out the details 
of the truck, as illustrated in Figure 43(b). Processing 
with the direct histogram specification method did not improve 
the image significantly, as can be seen in Figure 43(c). It 
should be kept in mind that this approach will never do any 
worse than histogram equalization, since histogram equaliza- 
tion is £ special case of the direct histogram specification 
method. The pertinent histograms are in Figures 43(d) 
through 43(g). 




(a) 



Figure hi. Results of processing the truck in Figure A0(a) 
with the histogram shown in Figure 41. 
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(a) 


(b) 



(c) 


Figure 43. Result of processing a low contrast picture of a 
truck . 
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Result No. 6 

As a final result, consider Figure AA . Figure AA(a) 
is the image of a woman's face covered with a shadow. This 
picture was taken from a photograph placed in a shadow. 

Note that the histogram equalization of this image produces 
a very noisy result, as seen in Figure AA(b) . The direct 
histogram specification method was effective in decreasing 
the prominence of the noise, as is illustrated in Figure 
A4(c) . The relevant histograms for this example are in 
Figures AA(d) through AA(g). 
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Figure 44. Result of processing the image of a woman’s 
face . 



CHAPTER V 


CONCLUSIONS 

The direct histogram specification method developed in 
Chapters II and III has been shovjn to be a valuable tool for 
interactive image enhancement. As was indicated in Chapter 
IV, this method can often improve the results obtainable 
with the v;ell“known hi.stogram equalization approach. Some 
experiments were described v.’here histogram equalization yielded 
enhanced images that were quite acceptable. Since this 
method is a special case of the al.goritlmi developed in this 
paper, however, the direct specification method will alvjays 
do at least as well as the classical equalization approach. 

The technique described in Chapter III for interactively 
specifying a histogram is just one of many that could have 
been implemented. Although the results obtained vjith this 
approach have proved satisfactory, other methods should be 
investigated and compared in terms of their interactive fea- 
tures. If this algorithm is to be implemented in a real-time 
system, the specification method needs also to be investigated 
with respect to its computational complexity. Its present 
form is suitable for the type of processing done here, but it 
may need to be simplified somewhat, in order to decrease the 
computation time required. 
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The histogram smoothing preprocessor presented in 
Chapter II has proven useful in reducing the error between 
the specified and resulting image histograms. Experiments 
have indicated, hov;ever, that the preproces.sed images tliat 
result from the histogram specification enhancement do not 
differ radically from those tliat were not preproces.sed. The 
preprocessor does not appear to have much impact on the 
final visual results. Theoretically, ho'wcver , it i.s of great 
significance if a uniform histogram is desired. Further work 
in this area could include the possibility of weighting the 
points in the neighborhood v.’hen performing the summation in 
Equation (16), page 

Ultimately, thi.s algorithm could be implemented in :\ 
real-time or almost real-time interactive image enhancement 
system. There are four p];incipal areas that must be 
considered : 

(1) calculation of the image hi.stogram, 

(2) computation of the transformation function, 

(3) computation of the new pixel values, and 

(4) reconstruction of tlic image for display. 

There are two approaches that can be taken. One is to use a 
digital system, the other a hybrid system. 

To proce.ss the image digitally (like all of the experi- 
ments presented in this paper) , the first three areas men- 
tioned above can be implemented with a special purpose 
computer. The fourth area involves the specification of a 
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system that will convert the digital data to an analog video 
signal for display on a monitor. The major problem with 
this approach would probably be the time required to retrieve 
the original image from storage and output the new image to 
the display device. This lias been the limiting faceor in 
the v7ork reported in tliis paper v/hen attemptinj' to .simulate 
a real-time system. 


The most attractive approach is to proce.ss the video 
signal in a hybrid system. This alleviates the problem oi 
reconstructing the imag,e , as v,'cl 1 as the cornputat ion of the 
new pixel values. The idea would be to generate an analog 
transfer function that is digitally controlled, as illustra- 
ted in Figure 45. The transfer function could be made up 
of a network of diode function generators that generate a 
piece-v.’ise linear approximation to the transformation 
function . 

A diode function generator has an i.nput-output voltag,e 
characteristic as illusti'ated in 1 if’.ure A6. ihe breakpoint, 
Eb, is variable as vjell as the slope, S. In a tran.sfer func- 
tion generator, the breakpoints could be fixed based on the 
number of specified seg^ments. Then the slope of each seg,- 
ment could be controlled digitally, using either a servo-con- 
trolled potentiometer or a solid state digital coefficient 


unit . 

The digital control could easily be supplied by an 
off-the-shelf microproco.ssor . Its function would be to 













Figure 46. Input-output characteristics of a diode function 
generator , 
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sample a part of the video sif^nal to compute a histogram 
representative of the image at that point in time. Using 
this information, a transforma tion function can then be 
computed and the .slope coefficients in the transfer function 
can be set. 

The hybrid approach has several advantages. First, it 
does not require a complex display system to generate a video 
signal, since it operates on the video directly. Second, 
since the picture does not need to be stored, the rend/vnite 
memory requirements for the system are minimal. The algorirlrn 
can be stored in read-only memory. The only random ac:cess 
memory required v;ould be to store the histogram and the coai- 
puted values of the transfer function. This requirement 
could be simplified considerably by dccrc-asing the spatial 
resolution and nuTnber of gray levels when obtaining the 
image histogram. 

In sutTUTiary, an image enhancement olgorit.hm has been 
developed and implemented. This interactive algorithm, 
which is based on directly specifying a histogram, has been 
shown to yield results that are often superior to those 
obtainable by the popular histogram oqualixation tccliniqiie. 

An algorithm has also been implemonted for improving the 
uniformity of the histogram. Finally, two methods have been 
proposed for implementing the direct histogram specif icatiion 
method of enhancement in a real-time interactive image 
enhancement system. 
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INTERACTIVE IMAGE ENIIANCEMENT SYSTEM PROGRAM 

This appendix describes the interactive image enhance- 
ment system program implemented in conjunction with this 

paper. Figure 47 illustrates a simplified flow chart of the 
program logic. 

Ihis program was written to run under the RT-11 Operat- 
ing System on the D.E.C. PDP 11/40. The main program and 
most of the subroutines were written in RT-11 FORTRAN IV. 

The other subroutines were written in PDP 11/40 assembler 
language . 

Three system subroutines were used. 'ASSIGN’ is used 
in conjunction with the DEFINE FILE statement to open a cer- 
tain direct access file. 'DATE' returns the current date in 
as ASCII character string. ' IDATE ' returns the current date 
in three integers denoting month, day and year. 

The program listing follows Figure 47. The main program 
is listed first followed by the FORTRAN subroutines in alpha- 
betical order. Finally, the two assembler language subroutine 
are listed. 
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START J 



Figure i\l . Simpl ified flov; chart nf interactive image 

enhancement aysteni pragram. 



COMPUTE 

HISTOGR/J-1 



Figure 47. (continued) 






Figure A7. (continued) 









Figure 47 . 


(continued) 







Figure 47 . 


(continued) 
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THIS PROGRAM IS INTENDED TO PERFORM AS AN 
INTERACTIVE IMAGE ENHANCEMENT SYSTEM (IIES) , 

CAPABLE OF PROVIDING THE FOLLOWING SERVICES: 

1) HISTOGRAM EQUALIZATION, 

2) HISTOGRAM FLATTENING , 

5) HISTOGRAM SPECIFICATION ENHANCEMENT, 

4) UTILITIES, 

IT PERFORMS IN A VERY SIMPLE QUESTION AND ANSWER 
FORMAT WITH THE USER CAPABLE OF SPECIFYING THE 
PRINTING OF HISTOGRAMS AND/OR STATISTICS, PLOTTING 
OF HISTOGRAMS , PLOTTING TR^J^S FORMAT IONS , ETC . 

OTHER ALGORITHMS CAN BE EASILY ADDED TO IT SINCE 
THE MAIN PROGRAM CONSISTS OF MOSTLY SUBROUTINE 
CALLS THAT PERFORM THE VARIOUS BASIC FUNCTIONS. 

EACH PROGRAM SECTION CONSISTS OF A SECTION OF 
OPERATOR DIALOG THAT ALLOWS THE USER TO SELECT 
CERTAIN OPTIONS, AS WELL AS THE CODE THAT ACTUALLY 
PERFORMS THE FUNCTION. THIS PROGRAM IS SET UP TO 
INTERCEPT ALL DIRECT ACCESS I/O ERRORS. IT PRINTS 
A MESSAGE ON THE CONSOLE DEVICE AND MAINTAINS 
CONTROL WITHIN THE PROGRAM. 

SUBROUTINES PLOTS, NEWPEN AND VIDPLT ARE THOSE THAT 
PERFORM PLOTTING ON THE VIDEO DISK. THEY ARE 
CALLED IN THE FOLLOWING MANNER: 

CALL PLOTS (IBK) ; WHERE IBK IS ON THE INTERVAL 
(0 , 15) , 0 REPRESENTING BLACK AND 15 
REPRESENTING WHITE. THIS SUBROUTINE ERASES 
THE SCREEN TO THE COLOR SPECIFIED BY IBK, AND 
INITIALIZES THE PEN TO (0,0). THE PEN COLOR 
IS SET TO VmiTE. 

CALL NEWPEN (IPEN) , WHERE IPEN IS ON THE 
RANGE (0,15). THIS CHANGES THE PEN COLOR. 

GALL VIDPLT (IX,IY,IP); WliERE (IX, lY) ARE THE 
X-Y COORDINATES WHERE THE PEN IS TO BE MOVED 
TO. IF IP IS +2, THE PEN WILL BE DOWN. IF IT 
IS +3, IT WILL BE UP. IF IP IS -2, THE PEN 
WILL BE DOWN AND A NEW ORIGIN DEFINED AT THAT 
POINT. IF IP IS -3, THE PEN WILL BE UP AND 
A NEW ORIGIN DEI’INED. 

SUBROUTINES USED - CAMERA, CRVPLT,GM1HX,GRAYEQ, HIST, 

HSTCRV , HSTPLT , HSTPRT . HSTSTF , 
INVCRV , POTENT . RECNUM . R. ^^ERR , 
RWLBLO , STAT , TRHIST , UNIRI^i.S , 
WRITEO,XLATE 
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THIS SECTION DEFINES THE ARRAYS THAT ARE USED 
THROUGHOUT THE PROGRMI. TO CONSERVE MEMORY, AS 
MANY ARRAYS AS ARE NEEDED ARE PASSED TO THE 
SUBROUTINES AS ARGUMENTS, EVEN IF THEY ARE NOT USED 
AS INPUTS TO OR OUTPUTS FROM THi: SUBROUTINE. THIS 
SECTION ALSO PERFORMS THE DIALOG NECESSARY TO 
DETERMINE WHICH ROUTINE IS DESIRED, AND TRANSFERS 
CONTROL TO THE APPROPRIATE PLACE. 


INTEGER*2 SYSDAT (1) . IPIC (2048) , PLTNUM,HSTPAS . CDFPAS 
INTEGER*2 ERFLG,XTAB(256) ,YTAB(256) ,ETAB(256) 

REAL*4 XH1ST(256) ,XT(11) ,YHIST(256) ,YT(U) ,ZHIST(256) 
REAL*4 ZT(ll) ,iHlST(256) ,TT(11) ,ST(4) ,GRHIST(230A) 
LOGICAL*! YES,NO,DAT(9) ,ANS,INFIL(14) ,OUTFIL(14) 
LOGICAL*! LPIC(4096) ,NAM!(8) 

LOGICAL TRUE , FALSE . PRS , PRH . PLH , PLC , TLD , GNC , POTS , OUTPUT 
LOGICAL OLDPRH 

DATA YES, NO, TRUE, FALSE/ 'Y' , 'N' . .TRUE. , .FALSE./ 
EQUIVALENCE (IPIC(l) ,LPIC(1)) . (GRHIST(!) , I?IC(!537) ) , 
1(GRHIST(1025) .YHIST(l)) ,(GRHIST(1281).ZHIST(!)) , 
2(GRHIST(l'i37) ,THIST(1)) , (GRHIST(!793) ,XTAB(1)) , 
3(GRHIST(1921) ,YTAB(!)) 

CALL DATE (DAT) 

CALL IDATE (I,J,K) 

S YSDAT ( ! ) « ( K - 7 2 ) +3 2* ( J+3 2*1 ) 

OLDPRH-FALSE 
WRITE (7,!) 

100 WRITE (7,2) 

READ (5,3) ANS 

IF (ANS. EQ. YES) GO TO 400 

WRITE (7,4) 

READ (5,3) ANS 

IF (ANS. EQ. YES) GO TO 700 

WRITE (7,5) 

READ (5,3) ANS 

IF (ANS. EQ. YES) GO TO 1000 

WRITE (7,6) 

READ (5,3) ANS 

IF (ANS. EQ. YES) GO TO 200 

STOP 'IIES FINISHED' 


THIS SECTION OF CODE IS THE UTILITY PROCESSOR. IT 
CAN PLOT AND PRINT HISTOGRAMS, PRINT STATISTICS 
AND PLOT CUMULATIVE DISTRIBUTION FUNCTIONS. 


200 WRITE (7.7) 
HSTPAS=0 
PLTNUM“0 
2!0 WRITE (7,8) 



i7 

READ (5.3) INFIX* 

IF CINFIL(l) .lQ.3a) GO TO 9999 
CAI*L ASaiGH (X , INFIL , 14 , 'OLD * , * NC ’ , 1) 

E®FINE fII*E 1 (130.1024,0. II) 

CALL RECNUM (LPIC.NREC.NAMl. 1 .ERFLG) 

IF ( (NREC . EQ . 0) . OR. (ERFLG , M. 0) ) GO TO 300 
WRITE (7.9) 

READ (5,3) AM,' ■ , ' : 

FLH>*FALSE 

IF (ANS.E^.tES) FLH^TOJE 
IF (.N0T.PLH) GO TO 220 
WRITE' ■(7rX0)v':/ ■ 

READ (5,3) ANS 

IF (ANS. NE. YES) GO TO 220 

WRITE (7. II) 

READ (S,X2) PLTNUM 
FLTNUM-MINO (PLTNOM , 4) 

PLTNUM=MAXO (PLTNUM , 1 ) 

HSTI'AS«0 

220 HSTPAS-HSTPAS+1 
WRITE (7,13) 

READ (5, 3) ANS 
PRR«FALSE 

IF (ANS.EQ.YES) PRH=TRUE 

17 (( .NOT.OLDFRH) .AND .PRS) WRITE (6 . 16) 

WRITE (7,14) 

READ (5,3) ANS 
PRS«FALSE 

IF (ANS.EQ.YES) PRS=TRtm 
WRITE (7.21) 

READ (5,3) ANS 
PLC«FALSE 

IF (ANS.EQ.YES) PLC=TRUE 
IF (.NOT.PLC) GO TO 230 
WRITE (7,10) 

READ (5.3) ANS 
CDFPAS=0 

IF (ANS.EQ.YES) GDFPAS«l 
230 CALL HIST (NREC,1,ERFU3,XPIG,XHIST) 

IF (ERFLG.NE.O) GO TO 300 
CALL STAT (XT.XHIST) 

CALL UNIRMS (NREC .XHIST.ERR) 

240 IF (.NOT. (PRS.OR.PEH)) GO TO 260 
250 WRITE (6.17) MMl, DAT 
260 IF (.NOT.PLH) GO TO 270 

CALL HSTPLT (PLTNBK.SSTPAS.XHIST) 

270 IF (.NOT.FRS) GO TO 280 
miTE (6,I8> XT 
WRITE (6,25) ERR 
WRITE ■ (6.X5) •: . ' . 

280 IF (.NOT.PRH) O) TO 290 


GALL ^TPRT (XHIST) 
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WRITE (6,16) 

290 IF (.NOT.PLC) GO TO 300 

CALL HSTCRV (NREC .XTAB .XlIIST) 

IF (PLH) PAUSE 'TYPE C/R FOR PLOT OF C.D.F.' 
CALL CRVPLT (CDFPAS , 2 ,XTAii) 

300 ENDFILE 1 
OLDPRH=PRH 
WRITE (7,19) 

READ (5,3) ANS 

IF (ANS. EQ. YES) GO TO 100 

GO TO 210 


THIS SECTION IS THE HISTOGR/iM EQUALIZATION 
PROCESSOR. THE OPTIONS AVAILABLE ARE PLOTTING AND 
PRINTING or HIST0G1L\MS, PRINTING OF STATISTICS AMD 
PLOTTING OF THE CWjULATIVE DISIRIBUTION FUNCTION. 
FIRST THE LABEL RECORD IS COPIED, THEN THE 
HISTOGRAM COMPUTED. THEN THE CUMULATIVE 
DISTRIBUTION FUNCTION IS CALCULATED AND THE PICTURE 
IS TRANSLATED. MEAhT-JHILE , THE OPTIONS ARE 
PERFORMi*:d AS requested. then THE NEW HISTOGRAM IS 
STORED IN THE LABEL RECORD. 


AOO WRITE (7,20) 

410 WRITE (7.9) 

•READ (5,3) ANS 
PT.H=FAT t’.P 

IF (ANS. EQ. YES) PLH=TRUE 
T^RITE (7,13) 

READ (5,3) ANS 
PRH-FALSE 

IF (ANS. EQ. YES) PRH=TRUE 
IF ((.NOT.OLDPRH) .AND.PRH) miTE (6,16) 
imiTE (7,14) 

READ (5,3) ANS 
PRS=FALSE 

IF (ANS. EQ. YES) PRS=TRUE 
\^RITE (7,21) 

READ (5,3) ANS 
PLC=FALSE 

IF (ANS.EQ.YES) PLC=TRUE 
420 WRITE (7,8) 

READ (5,3) INFIL 
IF (INFIL(l) .EQ.32) GO TO 9999 
CALL ASSIGN (1 , INFIL , 14 . 'OLD NC ’, 1) 
DEFINE FILE 1 (130, 1024 ,U, II) 

CALL RF.CNUM (LPIC ,NREC , N/\M1 , 1 , ERFLG) 

IF ((NREC. EQ.O) .OR. (ERFLG. NE.O)) GO TO 440 
WRITE (7,22) 

READ (5,3) OUTFIL 



IF (OUTFIL(l) .E0.32) GO TO 440 

CALL ASSIGN (2,dUTFIL,14, 'NEW , 'NCM) 

DEFINE FILE 2 (KRKC , 1024 . U , 12) 

CALL RWLBLO (1,2,’$', SYSDAT , LF I C , ERFLG) 

IF (ERFLG. NE.O) GO TO ^30 

CALL HIST (NREC.l, ERFLG, IPIC.XIIIST) 

IF (ERFLG. EQ.O) GO TO 4^0 

430 ENDFILE 2 

440 ENDFILE 1 
GO TO 550 

450 WRITE (7,23) 

CALL STAT (XT.XHIST) 

CALL HSTCRV (NREC ,XTAB ,XIIIST) 

CALL TRHIST (XTAB ,XH1ST , YKIST) 

\miTE (7,24) 

CALL STAT (YT,YHIST) 

CALL UNIRMS (NREC ,YHIST , ERR) 

IF (.NOT.PLH) GO TO 460 
CALL HSTPLT (2,1,XHIST) 

CALL HSTPLT (2,2,YHTST) 

460 CALL XIATE (NREC , 1 , 2 .ERFLG , IPIC ,XTAB) 

IF (ERFLG. NE.O) GO TO ^30 
CALL HSTSTF (2 ,YHIST .ERFLG . IPTC) 

IF (ERFLG. NE.O) GO TO 430 
ENDFILE 1 
ENDFILE 2 

470 IF (.NOT. (PRS.OR.PRH)) GO TO 540 

480 WRITE (6,17) NAM1.DAT 
WRITE (6,23) 

490 IF (.NOT.PRS) GO TO 500 
WRITE (6,18) XT 
WRITE (6,15) 

500 IF (.NOT.PRH) GO TO 510 
CALL HSTPRT (XHIST) 

WRITE (6,16) 

510 IF (.NOT. (PRS.OR.PRH)) GO TO 540 
WRITE (6,24) 

520 IF (.NOT.PRS) GO TO 530 
WRITE (6,18) YT 
WRITE (6,25) ERR 
WRITE (6,15) 

530 IF (.NOT.PRH) GO TO 540 
CALL HSTPRT (YHIST) 

WRITE (6,16) 

540 IF (.NOT.PLC) GO TO 550 

CALL HSTCRV (NREC ,YTAB , YHIST) 

IF (PLH) PAUSE ’TYPE C/R FOR PLOT OF C.D.F.’ 
CALL CRVPLT (1,2, XTAB) 

CALL CRVPLT (2,3,YTAB) 

550 WRITE (7,19) 

OLDPRH=PRH 
READ (5,3) ANS 
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IF (ANS.EQ.YES) GO TO 100 
WRITE (7,26) 

READ (5,3) ANS 
IF (ANS.EQ.YES) GO TO 410 
GO TO 420 


THIS SECTION OF CODE IS THE HISTOGRAM FL.ATTENING 
PREPROCESSOR. THE OPTIONS AVAILABLE ARE PRINTING 
AND PLOTTING 01’ HISTOGRAMS, PRINTING OF STATISTICS, 
AND THRESHOLDING. THE NUMBER OF ITERi\TIONS IS ALSO 
SELECTABLE. THE LABEL RECORD IS COPIED FIRST, THEN 
THE FIIES CLOSED AND REOPENED WITH THE LENGTH 
DEFINED TO BE TiiE LENGTH OF ONE ROW . THE 
SUBROUTINE 'GRAYEQ' IS CALLED AS MANY TIMES AS 
REQUESTED, AND THE OPTIONS ARE PERFORMED AS 
DESIRED. FINALLY, THE LAST HISTOGRAM IS STORED IN 
THE OUTPUT LABEL RECORD. 


700 WRITE (7,27) 

710 WRITE (7,9) 

READ (5,3) ANS 
PLH=FALSE 

IF (ANS.EQ.YES) PLH=TRUE 
WRITE (7,13) 

READ (5,3) ANS 
PRH=FALSE 

IF (ANS.EQ.YES) PRS-TRUE 

IF (( .NOT.OLDPRH) .AND.PRH) l^fRITE (6,16) 

WRITE (7,14) 

READ (5,3) ANS 

IF (ANS.EQ.YES) PRS-TRUE 
WRITE (7,38) 

READ (5,3) ANS 

IF (ANS.EQ.YES) TLD=TRUE 
WRITE (7,28) 

READ (5,12) ITER 
IF (ITER.LE.O) GO TO 900 
720 WRITE (7,8) 

READ (5,3) INFIL 
IF (INFIL(l) .EQ.32) GO TO 9999 
CALL ASSIGN (1 , INFIL , 14 ,' OLD' , ’NC ’, 1) 
DEFINE FILE 1 (130 , 1024 ,U , II) 

CALL RECNUM (LPIC ,NREC ,NAM1 , 1 , ERFLG) 

IF ((NREC.EQ.O) .OR. (ERFLG. NE.O)) GO TO 890 
WRITE (7,22) 

READ (5,3) OUTFIL 

IF (OUTFIL(l) .EQ.32) GO TO 890 

CALL ASSIGN (2 , OUTFIL, 14 , 'NEW' , 'NC ,1) 
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DEFINE FILE 2 (NREC , 1024 ,U , 12) 

CALL mrJLBLO (1 , 2 , ’ % ' , SYEDAT , LPIC . ERFLG) 

IF (ERFLG. NE.O) GO TO 880 

CALL HIST ( NREC , 1 , ERFLG .TIC, XHI ST) 

IF (ERFLG. NE.O) GO TO 880 
WRITE (7,23) 

CALL STAT (XT.XllIST) 

CALL UNIRMS (NREC ,XH 1ST , ERR) 

730 IF (.NOT. (PRS.OR.PRH) GO TO 750 
740 WRITE (6,17) NAM1.DAT 
750 IF (.NOT.PLH) GO TO 760 
CALL HSTPLT (4,1,XHIST) 

760 IF (.NOT.PRS) GO TO 770 
WRITE (6,23) 

WRITE (6,18) XT 
WRITE (6,25) ERR 
WRITE (6,15) 

770 IF (.NOT.PRH) GO TO 780 
CALL HSTPRT (XI i I ST) 

WRITE (6,16) 

780 CALL WRITEO (2 , NREC , IPIC , ERFLG) 

IF (ERFLG. NE.O) GO TO 880 
ENDFILE 1 
ENDFILE 2 

IF (NREC-34) 790,800,810 
790 NREC1=160 
NW=64 
GO TO 820 
800 NREC1=272 
NW=128 
GO TO 820 
810 NREC1=520 
NW=256 

820 DEFINE FILE 1 (NRECl , NW, U , II) 

CALL ASSIGN ( 2 , OUTFIL , 14 , ' OLD ’ , ’ NC ' , 1) 

DEFINE FILE 2 (NRECl , KEnLU . 12) 

DO 870 1=1, ITER 
ERFLG=I 

CALL GRAYEQ (GRHIST,NW,XHIST.TLD, ERFLG, IPIC, LPIC) 

IF (ERFLG. NE.O) TO GO 880 

WRITE (7,39) I 

CALL STAT (XT.XIIIST) 

CALL UNIRMS (NREC , XHIST, ERR) 

830 IF (.NOT.PLH) GO TO 840 
K=MOD(I,4)+l 
CALL HSTPLT (4, K, XHI ST) 

840 IF ( .NOT. (PRS.OR.PRH)) GO TO 870 
850 IF (.NOT.PRS) GO TO 860 
MITE (6,39) I 
WRITE (6,18) XT 
MITE (6,25) ERR 
WRITE (6,15) 
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860 IF (.NOT.PRH) GO TO 870 
CALL HSTPRT (XHIST) 

WRITE (6,16) 

870 CONTINUE 
ENDFILE 1 
ENDFILE 2 

CALL ASSIGN (2 , OUTFIL , 14 , 'OLD' , ’NC ,1) 
DEFINE FILE 2 (NREC , 102^1 ,U, II) 

CALL HSTSTF (2 .XllIST ,ERFLG , IPIC) 
ENDFILE 2 
GO TO 900 
880 ENDFILE 2 
890 ENDFILE 1 
900 WRITE (7,19) 

OLDPRH=PRlI 

READ (5.3) ANS 

IF (ANS. EO. YES) GO TO 100 

WRITE (7,26) 

READ (5,3) ANS 

IF (ANS. EQ. YES) GO TO 710 

GO TO 720 


THIS SECTION FUNCTIONS AS THE HISTOGRAM 
SPECIFICATION PROCESSOR. THE OPTIONS AVAILABLE ARE 
PRINTING AND PLOTTING OF HISTOGILAMS, PRINTING THE 
STATISTICS AND PLOTTING THE TR7vNS FORMATION 
FUNCTION. ALSO, IT IS POSSIBLE TO GENERATE ONLY A 
HISTOGRAM. THE USER HAS THE OPTION OF USING THE 
POTENTIOMETERS OR THE JOYSTICK-CURSOR AS INPUT. 
LABEL RECORD IS COPIED AND THE ORIGINAL HISTOGRAM 
AND C.D.F. COMPUTED. THE INVERSE IS TAKEN, THE 
FUNCTIONS COMBINED AND THE PICTURE IS TRANS FORl-lED . 
FINALLY THE NEW HISTOGRAM IS STORED IN THE LABEL 
RECORD. IF A NULL OUTPUT FILE NAME IS ENTERED, ALL 
OF THE ABOVE WILL HAPPEN EXCEPT THAT A NEW IMAGE IS 
NOT GENERATED. 


1000 WRITE (7,29) 

1010 miTE (7,30) 

READ (5,3) ANS 
GNC=FALSE 

IF (ANS. EQ. YES) GNC=TRUE 
WRITE (7.9) 

READ (5,3) ANS 
PLH=FALSE 

IF (ANS. EQ. YES) PLH=TRUE 
WRITE (7,13) 

READ (5,3) ANS 
PRH=FALSE 

IF (ANS. EQ. YES) PRH=TRUE 
IF ((.NOT.OLDPRH) .AND.PRH) 17RITE (6,16) 
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WRITE (7.1A) 

READ (5.3) ANS 
PRS- FALSE 

IF (ANS.EQ.YES) PRS-TRUE 
IF (GNC) GO TO 1240 
WRITE (7,37) 

READ (5.3) ANS 

IF (ANS*. EQ. YES) PLC=TRUE 
1020 WRITE (7,31) 

READ (5,3) ANS 
POTS-TRUE 

IF (ANS.EQ.YES) POTS-FAT, SE 
WRITE (7,8) 

READ (5,3) INFIL 

IF (INFIL (1 ) .EQ. 32) GO TO 9990 

CALL ASSIGN (1 , INFIL , 14 , 'OLD NC '. 1 ) 

DEFINE FILE 1 (130 , 1024 , U , II) 

CALL RECNUM (LPIC , NREC , NAIJI . 1 , EPFLG) 

IF ((NREC.EO.O) .OR. (ERFLG.NE.O)) GO TO 1230 
WRITE (7,22) 

READ (5,3) OUTFIL 
OUTPUT-TRUE 

IF (OUTFIL (1).EQ.32) OUTPUT-FALSE 
IF 0 ^30T . OUTPUT) GO TO 1030 
CALL ASSIGN ( 2 , OUTFIL , 14 , ' NEW ' , ' NC M) 
DEFINE FILE 2 (NREC , 1024 , U , 12) 

CALL RWLBLO (1 , 2 , ' A ’ , SYSDAT . LPIC , ERFLG) 

IF (ERFLG.NE.O) CO TO 1220 
1030 CALL NIST (NREC , i , ERFLG , IPIC .XHIST) 

IF (ERFLG.NE.O) GO TO 1220 
WRITE (7,23) 

CALL STAT (XT, XHIST) 

IF (.NOT.PLH) GO TO 1050 
1040 CALL HSTPLT (4,1, XHIST) 

1050 IF (.NOT. (PRS.OR.PRH)) GO TO 1080 
1060 VJRITE (6,17) NAM1.DAT 
TnFRITE (6,23) 

IF (.IsOT.PRS) GO TO 1070 
WRITE (6,18) XT 
WRITE (6,15) 

1070 IF (.NOT.PRH) GO TO 1080 
CALL HSTPRT (XHIST) 

WRITE (6,16) 

1080 CALL HSTCRV (NREC .XTAB , XHIST) 

CALL TRHIST (XTAB , XHIST ,THIST) 

WRITE (7,24) 

CALL STAT (TT,'rUIST) 

CALL UNIRMS (NREC , THIST, ERR) 

IF (.NOT.PLH) GO TO 1090 
CALL HSTPLT (4, 2, THIST) 



1090 IF ( .NOT. (PRB.OR.PRll)) GO TO 1110 
WRITE (6,2/<) 

IF (.NOT.PRS) GO TO 1110 
WRITE (6,18) TT 
WRITE (6,25) ERR 
WRITE (6,15) 

1100 IF (.NOT.PRH) GO TO 1110 
CALL IISTPRT (TUIST) 

WRITE (6,16) 

1110 IF (.NOT. POTS) GO TO 1120 
CALL POTENT (YIIIST , ST , NREC) 

GO TO 1130 

112.0 CALL CAMERA (YHIST , NREC , ZHIST , THIST) 

1130 WRITE (7,32) 

CALL STAT (YT. YHIST) 

IF (POT.S) WRITE (7,35) ST 
IF (.NOT.PLIl) GO TO llAO 
CALL HSTPLT (^i, 3. YHIST) 

1140 IF (.NOT. (PRS .OR.PRH)) GO TO 1160 
WRITE (6,32) 

IF (.NOT. PRS) GO TO 1150 
VJRITE (6,18) YT 
IF (.NOT. POTS) GO TO 1150 
WRITE (6,35) ST 

1150 WRITE (6,15) 

IF (.NOT.PRH) GO TO 1160 
CALL HSTPUT (YHIST) 

WRITE (6,16) 

1160 CALL HSTCRV (NRIX , YTAB , YHl ST) 

CALL IN VCRV (YTAL.ZTAB) 

CALL GMIHX (XTAB , YTAB , ZTAR) 

CALL TRHIST (2TAB , XHIST , ZPIST) 

WRITE (7,33) 

CALL STAT (ZT. ZHIST) 

CALL RMRERR (NREC , YHIST , ZHIST , ERR) 

IF (.NOT.PLH) GO TO 1170 
CALL HSTPLT (4, 4, ZHIST) 

1170 IF (.NOT. (PRS .OR.PRH)) GO TO 1190 
WRITE (6,33) 

IF (.NOT, PRS) GO TO 1180 
WRITE (6,18) ZT 
WRITE (6.3A) ERR 
WRITE (6,15) 

1180 IF (.NOT.PRH) GO TO 3190 
CALL HSTPRT (ZHIST) 

WRITE (6,16) 

1190 IF (.NOT. OUTPUT) GO TO 1200 
WRITE (7,36) 

READ (5,3) ANS 

IF (ANS. NE. YES) GO TO 1030 

CALL XLATE (NREC , 1 , 2 , ERFLG , IPIC , ZTAB) 

IF ( ERFLG. NE.O) GO TO 1220 
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CALL HSTSTF (2 , ZHIST, ERFLG , IPIC) 

ENDFILE 2 
1200 ENDFILE 1 

IF (.NOT.PLC) GO TO 1210 

IF (PLH) PAUSE ’TYPE C/R FOR PLOT OF CURVE’ 
CALL CRVPLT (1,3,ZTAB) 

1210 WRITE (7,19) 

OLDPRH-PRH 

READ (5,3) ANS 

IF (ANS. EQ. YES) GO TO 100 

WRITE (7,26) 

READ (5,3) ANS 

IF (ANS. EQ. YES) GO TO 1010 

GO TO 1020 

1220 IF (OUTPUT) ENDFILE 2 
1230 ENDFILE 1 
GO TO 1210 
1240 K-0 

1250 WRITE (7,31) 

READ (5,3) ANS 
POTS=TRUE 

IF (ANS. EQ. YES) POTS^FALSE 
IF (.NOT. POTS) GO TO 1260 
CALL POTENT (XllIST , ST , 10) 

GO TO 1270 

1260 CALL CAJffiRA (XHIST. 10 , ZHIST, THIST) 

1270 CALL STAT (XT, XHIST) 

IF (POTS) WRITE (7,35) ST 
IF (.NOT. PLH) GO TO 1290 
K-M0D(K,4)+1 

1280 CALL HSTPLT (4, K, XHIST) 

1290 IF ( .NOT. (PRS.OR.PRIl)) GO TO 1320 
WRITE (6,32) 

1300 IF (.NOT.PRS) GO TO 1310 
WRITE (6,18) XT 
IF (POTS) VmiTE (6,35) ST 
WRITE (6,15) 

1310 IF (.NOT.PRH) GO TO 1320 
CALL HSTPRT (XHIST) 

WRITE (6.16) 

1320 WRITE (7,19) 

OLDPRH=PRH 

READ (5,3) ANS 

IF (ANS. EQ. YES) GO TO 100 

WRITE (7,26) 

READ (5,3) ANS 

IF (ANS. EQ. YES) GO TO 1010 

GO TO 1250 
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1 

2 

3 

A 

5 

6 

7 

8 
9 

10 

11 

12 

13 

lA 

15 

16 

17 

18 


FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
1' MA = 


(’ INTERACTIVE IMAGE ENHANCEMENT SYSTEM*,//) 
c HISTOGRAM EQUALIZATION? ’, $) 

(lAAl) 

( ’ HISTOGRAM PREPROCESSOR? * , $) 

(' HISTOGRAM SPECIFICATION? ’, $) 

(’ UTILITY? ’,$) 

('OPICTURE UTILITY PROCESSOR’) 

(’ TYPE INPUT FILE NAME -- ',$) 

( ' PLOT HISTOGRAM? ’ , $) 

(• TIO,' NEW PLOT?', $) ^ 

(’ ', TIO, 'HOW MANY PLOTS ON SCREEN (1.2.A)?',$) 
(110) 

( ' PRINT HISTOGRAM? ’ , $) 

(’ PRINT STATISTICS? ' ,$) 

(IHO) 

r^PICTURE NAIIE -- ' , 8A1 , 20X, 'DATE: ',9A1) 

(' Ml = ',E11.A,' M2 = ’.Ell.A,’ M3 = ’,E11.A, 
',E11.A,/,' MODE “ ',F8.1,' AT ',F5.1. 


2' MEDIAN = ’,F5.1,' MIN = ’,F5.1,' MAX = ',F5.1, 


3' A3 = 
19 FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FOimT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORI^AT 
1 ’ H = 
36 FORJIAT 
FORMAT 
FORMAT 
FORMAT 
STOP 
END 


,E11.A) 

’ .$) 


20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 


37 

38 

39 
9999 


Ell. 4 ' A4 = 

(' SwiTCH MODES? . , , 

rOHISTOGRAM EQUALIZATION PROCESSOR’) 


PLOT CUl-IULATIVE DISTRIBUTION FUNCTION? ’, $) 
TYPE OUTPUT FILE NAME -- ’,$) 

ORIGINAL IMAGE’) 

HISTOGRAM EQUALIZED IMAGE') 

RMS ERROR WITH RESPECT TO UNIFORM = ’,E11.4) 
CHANGE PARAMETERS? ' ,$) 

OHISTOGRAM FLATTENING PREPROCESSOR’) 

ITERATION COUNT? ' , $) 

ODIRECT HISTOGRAM SPECIFICATION PROCESSOR ) 
GENERATE HISTOGRAM ONLY?',$) 

INPUT FROM JOYSTICK-CURSOR? ' ,$) 

DESIRED IMAGE’) 

RESULTING IMAGE’) 

RMS ERROR WITH RESPECT TO DESIRED = ’,E11.4) 
M = ',E11.4,’ LS = '.Ell. 4,’ RS = ’,E11.4. 
,E11.4) 

( ’ SATISFIED? ' , $) 

(’ PLOT TRANS F6R^tATION CURVE? ’,$) 

(' THRESHOLD? ' .$) 

(' ITERATION ^-'',14) 

IIES FINISHED' 


( 

(' 

(’ 

(' 

(’ 

( 

( 

(’ 

( 

(■ 

(' 

(* 

(’ 

(' 

(’ 
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SUBROUTINE CAMERA (HISTO.NREC ,XT, YT) 


THIS SUBROUTINE GENERATES A HISTOGRAM ’HISTO* THAT 
IS ENTERED VIA THE DIG ITIEING CAMERA USING THE 
JOYSTICK-'amSOR TO TRACE THE CURVE. THE HISTOGRAM 
IS SCALED TO THE CORRECT NUMBER OF POINTS USING 
•NREC’, THE NUMBER OF 2048 BYTE RECORDS IN THE 
FILE. ’XT' AND ’YT' ARE USED AS WORK SPACE. 


THE USER FIRST POSITIONS THE CURSOR ALONG THE ZERO 
LINE OF THE HISTOGPw\M. THIS VALUE IS ENTERED AS 
THE ZERO Y VALUE. THEN THE MAXIMITI X VALUE IS 
ESTABLISHED. THE USER HAS THE OPTION OF DIGITIZING 
CONTINUOUSLY OR PIECE-WISE. IF PIECE -WISE. EACH 
POINT IS SAMPLED BY TYPINC? A CARKIAGi: RKTURM . THE 
X-Y COORDINATES ARE SAMPLED AND S'lOKED IN 'XT’ AND 
'YT' RESPECTIVELY. ANY POINT l^fHERE Y IS LESS THAN 
THE ZERO REFERENCE ARE IGNORED. ANY POINTS WIERE 
X IS LESS THAN A PREVIOUS VALUE ARi: ALSO DISCARDED. 
WHEN X REACHES XMAX, THE DIGITIZATION IS COMPLETE. 
THE X VALUES ARE THEN RESCALED TO (1,256). THE 
Y VALUES ARE THEN USED AS BREAKPCINTS TO GENERATE 
A PIECE-WISE LINEAR HISTOGRAM'. IF MORE THAN ONE 
ELEMENT IN X HAS THE S.A^IE VALUE. THE FIRST ONE IS 
USED. THE REST ARE DISCARDED. THE HISTOGRAM IS 
THEN RESCALED TO THE CORRECT NUMBER OF POINTS. 

SUBROUTINES USED - CURSOR 


1NTEGER*2 XT(512) ,YT(512) ,XZ,XM,YZ 
REAL^-4 HISTO(256) 

LOGICAL*! ANS.Y 
DATA Y/ 'Y'/ 

100 PAUSE 'ENTER ZERO Y' 

CALL CURSOR (IX, YZ) 

YZ“513 -YZ 

PAUSE 'ENTER MAXIMUM X' 

CALL CURSOR (XM,IY) 

WRITE (7,1) 

READ (5,2) ANS 

PAUSE 'POSITION CURSOR AT LEFT -MOST POINT AND TYPE CR' 
1=1 

IXT=-1 

110 IF (ANS. NE.Y) PAUSE 'TYPE C/R' 

CALL CURSOR (IX. lY) 

IF (IX.GT.XM) GO TO 120 
IY=511-IY 

IF ((lY.LT.YZ) .OR. (IX.LE.IXT)) GO TOllO 
IXT=IX 
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XT<I)»IX 

YT(I)»iy-YZ 

I«I+1 

IF (I.LE.512) GO TO 110 
120 I»I-1 

WRITE (7,3) I 
XM“XT(I) 

XZ«XT(1) 

DO 130 J»1,I 

XT(J)»(XT(J) -XZ)*255 . 0/ (XM-XZ)+1 . 5 
130 CONTINUE 
K»1 
J-1 
IXT==-1 
140 IX=XT(K) 

IF (IXT.EQ.IX) GO TO 170 
IF (IX.NE.J) GO TO 150 
HISTO(J)“YT(K) 

J=J+1 
GO TO 170 

150 FYD“YT(K)-HISTO(J-l) 

IXD=IX-J 

DO 160 M“1,IXD 

HISTO ( J ) -HI STO ( J - 1) +FYD/ IXD 

J-J+1 

160 CONTINUE 
170 IXT-IX 
K-K+1 

IF (K.LE.I) GO TO 140 
T-0. 0 

DO 180 I»l,256 
T=T+HIST0(I) 

180 CONTINUE 

FH»(NREC-2)*2048 . 0/T 
DO 190 1-1,256 
HISTO(I)»HISTO(I)*FH 
190 CONTINUE 

1 FOm-lAT (' CONTINUOUS OR PIECE-WISE*, 

1’ (Y FOR CONTINUOUS)?* ,$) 

2 FORMA? (Al) 

3 FORMAT r* DIGITIZATION COMPLETE IN *,I4,’ POINTS') 
RETURN 

END 


SUBROUTINE CRVGEN (HIST,IA,IB,IC.AA,BB,CC) 


THIS SUBROUTINE GENERATES THE PIECE-WISE LINEAR 
HISTOGRAM * HI ST*. *IA* , *1B* AND ’IC* ARE THE 
X COORDINATES OF THE THREE BREAKPOINTS. ’AA* , 
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C *BB* AND ‘CC’ ARE THE Y COORDINATES. THE X 

C COORDINATES ARE ON THE RANGE (1,256). THE 

C HISTOGRAM IS PIECE-WISE LINEAR ON THE FOLX^ING 

C SEGMENTSt (1,0) - (lA.AA) - (IB.BB) - (IC.CC) - 

C (256,0), IT COULD BE THAT ’IA'»1, OR »IC’«256. 

C IT CHECKS FOR THIS AND TAKES THE APPROPRIATE 

C ACTION. 

C 

C SUBROUTINES USED - NONE 

C 

c 

REALM HIST (2 56) 

IF (lA.EQ.l) GO TO 110 

DO 100 T^l.IA 

HIST (I)-(I-1)*AA/(IA-1) 

100 CONTINUE 
GO TO 120 
110 HIST(1)=AA 
120 K»IA'i-l 

IF (lA.NE.IB) GO TO 130 
HIST (IB)-BB 
GO TO 150 
130 DO 140 I=K,IB 

HIST (I)=liIST(IA) + (I+l-K)*(BB-AA)/(IB-IA) 

140 CONTINUE 
150 K»IB+1 

IF (IB.EQ.IC) GO TO 170 
DO 160 I=K,IC 

HIST(I)-HIST(IB)+(1+1-K)*(CC-BB)/(IG-IB) 

160 CONTINUE 

170 IF (1C.EQ.256) GO TO 190 
HIST(IC)=CC 
K-lC+1 

DO 180 I=K,256 

HIST(I)-HIST(lC)+(I+l-K)*(-CC)/(256-IC) 

180 CONTINUE 
190 RETURN 
END 



SUBROUTINE CRVPLT (M , REF . TABLE) 

C 

G 

C THIS SUBROUTINE PLOTS A TRANSFORMATION CURVE PASSED 

G TO IT IN 'TABLE'. IT ASSUMES THAT THE VALUES IN 

C 'TABLE' ARE IN THE RANGE (0,255). 'M' INDICATES 

C WHETHER OR NOT TO PLOT A GRID AND CLEAR SCREEN. 

C IF M IS 1, IT IS DONE. IF NOT, THE CURVE IS 

C PLOTTED OVER WHATEVER IS THERE. A 45 EEGREE 

C REFERENCE LINE IS PLOTTED IF 'REF '-2 AND ’M'»l 

C THIS SUBROUTINE IS SET UP SO THAT IF A HISTOGRAM IS 
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C PLOTTED USING ’HSTPLT' , WITH 'PASS‘^1, 

C THE C.D.F. CAN BE PLOTTED OVER IT WITH A CALL TO 

C 'CRVPLT' WITH ’M' NOT EQUAL TO 1. 

C SUBROUTINES USED - NEWPEN .PLOTS ,VIDPLT 


INTEGER*2 TABLE ( 2 5 6 ) , REF 
IF (M.NE.l) GO TO 120 
100 CALL PLOTS (0) 

CALL VIDPLT (17.13.-3) 
CALL NEWPEN (6) 

DO 110 1=1,9 
IY=50*(I-1) 

CALL VIDPLT (0,IY,3) 

CALL VIDPLT (512. lY, 2) 
IX=64*(I-1) 

CALL VIDPLT (IX. 0,3) 

CALL VIDPLT (IX. 4 00. 2) 
110 CONTINUE 

REF=IABS (REF) 

CALL VIDPLT (0.0, REF) 
CALL NEWPEN (12) 

120 CALL VIDPLT (0,0.3) 

DO 130 1=1,256 

IY=1 . 5625*TABLE(I)+0 .5 

IX=2’*f(I-l) 

CALL VIDPLT (IX.IY.2) 
130 CONTINUE 
RETURN 
END 


SUBROUTINE DIFF (DELT.IA.IB .IM.TLD.IC , J) 


THIS SUBROUTINE CO^IPUILS THE DIFFERENCE SUMl-tATION 
‘DELT’ IF ’IC*-1. 'lA' IS THE ORIGINAL POINT. 

»IB' IS THE 9 POINT SUMItATION TO REPLACE 'IA‘. 

IF 'TLD’ IS FALSE, THRESHOLDING IS NOT DESIRED SO 
IT RETURNS. IF ’IC'-2, THE DIFFERENCE IS COMPUTED 
AND TESTED TO SEE IF IT LIES WITHIN THE BOUNDARY 
AROUND 9*IA. IF NOT, 'J' IS USED TO REASSIGN THE 
POINT TO SOME VALUE WITHIN THE LIMITS. IM IS THE 
MEAN OF THE SUMMATION ’DELT’ IF ’IC’=2. 


SUBROUTINES USED - NONE 


REAL*4 DELT (2) 
LOGICAL TLD 
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IF (.NOT.TLD) RETUiy^ 

ID=9*IA 

I=IABS(ID-IB) 

IF (I.EQ.O) RETURN 
GO TO (100,110) ,IC 

100 DELT(1)=DELT(1)+I 
DELT(2)=DELT(2)+1.0 
RETURN 

110 IF (I.LE.IM) RETURN 
IS=ISIGN(1, (IB-ID)) 
J=MOD(J,IM+l)+l 
IB=ID+J*IS 
RETURN 
END 


SUBROUTINE EQUAL (IW , HIST , BELT , IM) 


THIS SUBROUTINE COMPUTES THE CUMULATIVE 
DISTRIBUTION FUNCTION OF THE HISTOGRAI4 ’HIST’ AND 
STORES IT BACK IN ’HIST’ WITH VALUES ON THE RANGE 
(0,255). IT ALSO COMPUTES THE MEAN OF THE 
DIFFERENCE SUMMATION 'DELT’ AND RETURNS IT IN IM . 
’NW’ IS THE NUMBER OF POINTS IN A ROW OF THE PICTURE. 


SUBROUTINES USED - NONE 


REAL*4 HIST(2304) ,DELT(2) 


IM=0 

WORD==FLOAT (W) *FI.OAT (Nl^) 

SCALE-255.0/WORD 

T^=0.0 


DO 100 1=1,2304 
T=T+HIST(I) 

HIST(I)=AINT(T*SCALE-l-0 . 5) -128 . 0 

IF (HIST (I). LT. 0.0) HIST(I)=HIST(I)+256.0 


100 CONTINUE 

IF (DELT (2) . EQ . 0 . 0) RETURN 
IM=DELT (1) /DELT (2)+0 . 5 
RETURN 
END 


SUBROUTINE FIX (IC,NREC,HIST,OPIC,L0PIC,NW,HISTl,ERFLG) 


IF ’IC’=1, THIS SUBROUTINE UPDATES THE HISTOGRAM 
•HIST' USING THE DATA IN 'OPIC. ’NW' IS THE 
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NUMBER OF DATA POINTS IN 'OPIC‘ . 'LOPIC' IS THE 
SAME ARRAY AS 'OPIC’ , BUT IN LOGICAL*! RATHER THAN 
INTEGER*2. IF ’IC'«2, THIS SUBROUTINE TRANSLATES 
THE DATA IN 'OPIC' USING ’HIST’ AS THE 
TRANSFORMATION. IT ALSO COMPUTES THE 256 POINT 
HISTOGRAM ’HIST!’. *OPIC IS THEN OUTPUT TO RECORD 
’NREC IN THE OUTPUT FILE. ’ERFLG’-O IF THERE WERE 
NO ERRORS . 

SUBROUTINES USED - NONE 


INTEGER*2 OPIC (512) ,ERFLG 
LOGICAL*! LOPIC (1024) 

REAL*4 HIST(2304) ,HIST1(256) 

ERFLG-0 

NB=NI^+NW 

GO TO (100,L20),IC 
100 DO 110 1=1, 

K=OPIC(I)+l 
HIST (K)=HIST(K) +1.0 
110 CONTINUE 
RETURN 

120 DO 130 1=1, N\^ 

K=OPIC(I)+l 

OPIC(I)=HIST(K) 

M=HIST(K) 

HISTl (M+1)=HIST1 (M+l)+l . 0 
130 CONTINUE 

WRITE (2'NREC,END=140,ERR=150) (LOPIC (I) ,1=1, NB, 2) 
RETURN • . / 

140 WRITE (7,1) 

GO TO 160 
150 WRITE (7,2) 

160 ERFLG=1 

1 FORMAT ( ' END OF FILE ON UNIT 2 IN ” FIX " ’ ) 

2 FORMAT (’ I/O ERROR ON UNIT 2 IN “FIX'") 

RETURN 

END 


SUBROUTINE GMIHX (H,G,T) 

C 

c 

r CREATES THE TRANSFORMATION 

O T«G(H(X)). 

c 

C SUBROUTINES USED - NONE 

C 

C 

INTEGER*2 H(256) ,G(256) ,T(256) 

DO 100 1=1,256 
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J=H(I)+1 
T(I)=G(J) 
100 CONTINUE 
RETURN 
END 


SUBROUTINE GR/\YEQ (KIST.NX.HISTl ,TLD,ERFLG, IPIC.LPIC) 


THIS SUBROUTINE PERFORMS THE HlSTOGRAI'l FLATTENING 
OPERATION. TWO COMPLETE PASSES ARE MADE THROUGH 
THE PICTURE. ON THE FIRST PASS, THE 2,304 POINT 
HISTOGRAM ’HIST' IS COMPUTED, AS WELL AS THE 
DIFFERENCE SUMMATION ' DIXT ' . THE C . D . F . OF ' HIST ' 
IS THEN COMPUTED AND STORED BACK INTO 'HIST'. THE 
SECOND PASS THEN USES 'HIST' TO TRANSLATE THE DATA 
AND THRESHOLD IT IF 'TLD' IS TRUE. 'HISTl' IS THE 
256 POINT HISTOGRAM OF THE RESULTING IMAGE. 
'ERI-’LG'=0 IF THERE WERE NO ERRORS. 'IPIC IS USED 
AS WORKSPACE. "LPIC IS THE SAME ARPvAY BUT IN 
LOGICAL*! FORM. 'NX' IS ONE KAT.F THE NUMBER OF 
POINTS IN ONE ROW. THE SUBROUTINE OPERATES ON 3 
ROWS AT A TIME TO COMPUTE A 9 POINT SmiMATION. 

THEY ARE STORED IN THE FIRST 3 ROWS OF ' IPIC ' . 

THE TOP AND BOTTOM ROWS AND ALONG THE SIDES ARE 
ONLY 6 POINT SUMl'iATIONS , BUT THEY ARE MULTIPLIED 
BY 1.5. 

THE FOUR CORNERS ARE 4 POINT SUMMATIONS MULTIPLIED 
BY 2.25. THE TOP AND BOTTOM ROWS ARE PROCESSED 
FIRST, THEN THE OTHERS. THE 1/0 IS SET UP SO THAT 
ONLY 1 ROW NEED BE READ AT A TIME, RATHER TIiAN ALL 
THREE. IT ASSUMES INPUT ON LUN 1, OUTPUT ON LUN 2. 

SUBROUTINES USED - DI7F, EQUAL, FIX, JUG 


INTEGER*2 IPIC (512, 4) ,OPIC(512) ,ERFLG 
LOGICAL*! LPIC (1024 , 4) , LOPIC (1024) 
REAL*4 HIST (2304) ,DELT(2) , HISTl (256) 
LOGICAL TLD 

EQUIVALENCE (OPIC(l) .LOPIC(l)) 

L=1 

IF (ERFLG.NE.l) L=2 

ERFLG-0 

NR=2048/NX 

NW=NX+NX 

NB=NW+NW 

IEN=NR+NW 

IEN1=^IEN-1 

NW1=NW-1 
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DO 100 I«l,2304 
HIST(I)=0,0 
100 CONTINUE 
DELT(1)=0.0 
DELT(2)=0.0 
DO 110 1=1,256 
HIST1(I)=0.0 
110 CONTINUE 

DO 120 J=l,3 
DO 120 1=1, 

IPIC(I,J)=0 
120 CONTINUE 

DO 180 IC=1,2 

READ (L*NR+1,END=190,ERR=200) (LPTC(J , 1) , J=1,NB, 2) 

READ (L’NR+2,END=190,ERR=200) (LPIC(J , 2) , J=1 ,NB, 2) 

CALL JUG (IPIC,NW,1) 

CALL JUG (IPIC,KW,2) 

0PIC(1)=(1PIC(1,1)+IPIC(1.2)+IPIC(2,1)+IPIC(2,2))*2.25 

IJK=0 

CALL DIFF (DELT,IPIC(1,1) ,OPIC(l) ,IM,TLD,IC,IJK) 
OPIC(mO*=(IPIC(NW,l)+lPIC(mi,l)+IPIC(NW,2) + 
1IPIC(NW1,2))*2.25 

CALL DIFF (DELT , IPIC (IW , 1) , OPIC (NW) , IM , TLD , IC , UK) 

DO 130 J=2,Nin 

OPIC ( J)= (IPIC (J-l , 1)+IPIC (J , D+IPIC (J+1 , D+IPIC ( J-1 , 2)4- 
IIPIC ( J , 2)+IPIC (J+1 , 2) ) *1 .5+0.5 
CALL DIFF (DELT,IPIC(J,1) ,OPIC(J) ,IM,TLD,IC,IJK 
130 CONTINUE 

CALL FIX (IC,NR+l,HIST,e^IC.LOPIC,NW,HISTl,ERFLG) 

IF (ERFLG.NE.O) RETURN 

READ (L’IEN1,END=190.ERR=200) (LPIC(J,1) , J=1,NB, 2) 

READ (L'IEN,END=190,ERR=200) (LPIC(J,2) , J=1,NB,2) 

CALL JUG (IPIC,NW,1) 

CALL JUG (IPIC,NW,2) 

OPIC(l)=(IPIC(l , 1)+IPIC(1 , 2)+IPIC(2 , 1)+IPIC(2 , 2) )*2. 25 
CALL DIFF (DELT,IPIC(1,2) ,OPIC(l) ,IM,TLD,IC,IJK) 

OPIC (NW) = ( IPIC (NW, D+IPIC (KWl , l)+IPIC(^Tv^ 2) + 

IIPIC (NT^l , 2)) *2. 25 

CALL DIFF (DELT , IPIC (NW, 2) , OPIC (NW) , IM, TLD , IC , UK) 

DO 140 J=2,m^Il 

OPIC(J) = (IPIC (J-1, D+IPIC (J,1)+IPIC(J+1, D+IPIC (J-1, 2) + 
IIPIC (J , 2)+IPIC(J+l , 2) )*1 .5+0.5 
CALL DIFF (DELT,IPIC(J,2) ,OPIC(J) ,IM,TLD,IC,IJK) 

140 CONTINUE 

CALL FIX (IC.IEN,HlST,OPIC,LOPIC,tR^,HISTl,ERFLG) 

IF (ERFLG.NE.O) RETURN 
K=NR+2 

READ (L‘K-1,END=190,ERR=200) (LPIC(J, 1) , J=1,NB, 2) 

READ (L'K,END=190,ERR=200) (LPIC(J,2) , J=1,NB,2) 

CALL JUG (IPIC,NW,D 
CALL JUG (IPIC,NW,2) 

IK=2 
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DO 170 M1=3,NW 
IR0W=M0D(IK,3)+1 

READ (L'K+1, END-1 90, ERR-200) (LPIC(J.IROW) , J-1,NB,2) 
lK-IK+1 

IROWI-MOD ( IROWH-1 , 3 ) +1 
CALL JUG (IPIC,NW,IROW) 

0PIC(1)=(1PIC(1,1)+IPIC(1,2)+IPIC(1,3)+IPIC(2.1)+ 

lIPIC(2,2)+IPIC(2,3))*1.5+0.5 

CALL DIFF (DELT,IPIC(l,IROWl) ,OPIC(l),IM,TLD,IC,IJK) 
one (NW) - (IPIC (NWl , 1) -MPIC (NWl , 2)+IPIC (NWl , 3 ) + 
IIPICCNW, 1)+IPIC(NW, 2)+IPIC(NW, 3^ )*1 .5+0.5 
CALL DIFF (DELT.IPIC(NW,1R0W1) ,0P1C(NW) ,IM,TLD,IC,IJK) 
DO 160 M2=2,NW1 
• OPIC(M2)-0 
DO 150 M3-I,3 
DO 150 M4=l,3 
K1-M2+M3-2 

OPIC (M2 ) -OPIC (M2) +IPIC (K1 , M4) 

150 CONTINUE 

CALL DIFF (DELT,IPIC(M2,IROWl),OPIC(M2).IM,TLD,IC.IJK) 
160 CONTINUE 

CALL FIX (IC,K, HIST, OPIC, L0PIC,NW,HIST1,ERFLG) 

IF (ERFLG.NE.0) RETURN 

170 CONTINUE 

IF (IC.EQ.2) GO TO 180 
CALL EQUAL (NW,HIST,DELT, IM) 

180 CONTINUE 
RETURN 

190 WRITE (7,1) 

GO TO 210 
200 WRITE (7,2) 

210 ERFLG-1 

1 FORMAT (' END OF FILE ON UNIT 1 IN ’’GRAYEQ'") 

2 FORMAT (’ I/O ERROR ON UNIT 1 IN "GRAYEQ"’) 

RETURN 

END 


SUBROUTINE HIST (NREC,LUN,ERFLG, IPIC,HISTO) 


THIS SUBROUTINE COMPUTES THE HISTOGRAM OF THE 
PICTURE IN FILE 'LUN' . THE HISTOGPvA.M IS STORED IN 
'HISTO* . IT FIRST CHECKS TO SEE IF THE HISTOGRAM 
HAS BEEN STORED IN THE LABEL RECORD. IF IT HAS, 

IT IS MOVED INTO ’HISTO' . IF IT HAS HOT, IT IS 
COMPUTED AND THEN \^RITTEN INTO THE LABEL RECORD, 
*IPIC IS USED AS AN I/O BUFFER. 'ERFLC' IS ZERO 
IF THERE WERE NO ERRORS. 'NREC IS THE NUl^ER OF 
RECORDS IN THE FILE. 


oooo 


SUBROUTINES USED - READO 


INTEGER*2 IPIC(2048) ,ERFLG,IHIST(512) 
REAL HISTOC256) ,TEMP(256) 

EQUIVALENCE (IHIST(l) ,TEMP(1)) 

ERELG=0 

READ (LUN’1,END«170,ERR-180) (IPIC(I),I= 
IF (IPIC(512) .NE.-l) GO TO 120 
DO 100 1-1,512 
K-512+I 

IHIST(I)=1PIC(K) 

100 CONTINUE 

DO 110 1=1,256 
HISTOCI)-TEMP(I) 

110 CONTINUE 
RETURN 

120 DO 130 1=1,256 
HISTO(I)=0.0 
130 CONTINUE 

DO 140 I=3,NREC 

CALL READO (LUN,I,ERFLG,IPIC) 

IF (ERFLG.NE.O) GO TO 999 
DO 140 J=l,2048 
M=IP1C(J)+1 
HISTO (M) =HISTO (M) +1 . 0 
140 CONTINUE 

READ (LUN'1,END=170,EKR=180) (IPIC(I),I= 
IPIC(512)=-1 
DO 150 1=1,256 
TEMP (I) -HISTO (I) 

150 CONTINUE 

DO 160 1=1,512 

V— CT 0-LT 

IPIC(K)=IHIST(I) 

160 CONTINUE 

WRITE (LUN’1,END=170,ERR=180) (IPIC(I), 
RETURN 

170 miTE (7,1) LUN 
GO TO 190 

180 WRITE (7,2) LUN 
190 ERFLG-1 

1 FORMAT (’ END OF FILE ON UNIT ‘,13,’ IN 

2 FOI^tAT (' I/O ERROR ON UNIT ‘,13,’ IN ’ 
999 RETURN 

• END 


1,1024) 


1,1024) 


=1,1024) 


• ‘HIST' ") 
'HIST' ' *) 
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SUBROUTINE HSTCRV (NREC, TABLE, HISTO) 


THIS SUBROUTINE COMPUTES THE CUMULATIVE 
DISTRIBUTION FUNCTION 'TABLE' OF THE HISTOGRAM 
'HISTO'. IT IS SCALED ON THE RANGE (0.2b')) SO THxVf 
IT MAY RE USED AS A TRANSFOl^IATION TABLE FOR 
HTSTOGRAII EQUALIZATION. 'NREC' DENOTES THE NUMBER 
OF 20A8 BYTE RECORDS IN THE FILE. 

SUBROUTINES USED - NONE 


INTEGER*2 TABLE (256) 
REAL HISTO (256) 
PlXEL^-2048 . 0*(NREC-2) 
SCALE=255.0/ PIXEL 
TEMP=0.0 
DO 100 1=1,255 
TEMP=SCALE*1IIST0 (I ) +TEMP 
TABLE(I)“TEKP+0,5 
100 CONTINUE 

TABLE (256) =255 

RETURN 

END 


SUBROUTINE HSTPLT (NUM, PAS S , HISTO) 


THIS SUBROUTINE PLOTS A SERIES OF HISTOGR/lMS ON THE 
TELEVISION. 'NUM' INDICATES THE NUMBER OF 
HISTOGRAMS TO BE PLOTTED ON THE SCREEN AT ONCE - 
EITHER 1, 2, OR 4. 'PASS' INDICATES WHICH 
HISTOGRAM IS CURRENTLY BEING PLOTTED. 'HISTO' IS 
THE HIETOGR.\M TO BE PLOTTED. IF 'PASS'=1, THE 
SCREEN IS CLEARED AND THE GRID PLOTTED. SUBSEQUENT 
PASSES MERELY HOVE THE ORIGIN UP THE LEFT HAND 
VERTICAL AXIS. 

SUBROUTINES USED - NEWPEN .PLOTS .VIDPLT 


INTEGER*2 NUM. PASS 
REAL*4 HISTO (2 5 6) 

IF (NUM.EQ.3) NUM=4 
INC=400/NUM 

IF (PASS.NE.l) GO TO 100 
CALL PLOTS (0) 

100 CALL NEWPEN (6) 

IF (PASS.NE.l) GO TO 120 
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CALL VIDPLT (17, 3, -3) 

DO 110 1=1,5 
IY=100*(I-1) 

CALL VIDPLT (0,TY,3) 

CALL VIDPLT (515, I Y, 2) 

CALL VIDPLT (515,IY+1,2) 

CALL VIDPLT (0,IY+1,2) 

110 CONTINUE 

CALL VIDPLT (0,0,2) 

CALL VIDPLT (1,0,2) 

CALL VIDPLT (1,401,2) 

CALL VIDPLT (514,0,3) 

CALL VIDPLT (514,401,2) 

CALL VIDPLT (515,401,2) 

CALL VIDPLT (515,0,2) 

CALL VIDPLT (2,1, -3) 

120 IF (PASS.EQ.l) GO TO 130 
CALL VIDPLT (0,INC,-3) 

130 DO 140 1=1,3 
IX=128*I 

CALL VIDPLT (IX, 3, 3) 

CALL VIDPLT (IX, -3, 2) 

CALL VIDPLT (IX+1,-3,2) 

CALL VIDPLT (IX+1,3,2) 

140 CONTINUE 

CALL VIDPLT (0,0,3) 

FM=0 . 0 

DO 150 1=1,256 

IF (FM.LT.H1ST0(D) FM=HISTO(I) 
150 CONTINUE 

SCALE=(TNC-1)/FM 
CALL NEWPEN (9) 

DO 160 1=1,256 

IF (HISTO(I) .EQ.0.0) GO TO 160 
IX=2*(I-i) 

CALL VIDPLT (IX. 0,3) 

IY=HISTO (I) *SCALE+0 . 5 
CALL VIDPLT (IX,IY,2) 

160 CONTINUE 
RETURN 
END 


SUBROUTINE HSTPRT (HIST) 

THIS SUBROUTINE TAKES A 256 POINT HISTOGRAM 'HIST' 
AND PLOTS IT ON THE LINE PRINTER. 

SUBROUTINES USED - NONE 
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REAL*4 HIST (256) 

LOGICAL*! LINE (101) 

DATA LlNE/101*’*’/ 

T«0 . 0 

DO 100 1=1,256 
T=AMAX1(T,HIST(I)) 

100 CONTINUE 
T=100.0/T 
DO 110 1=1,256 
J=T*HIST(I)+1.5 
K=I-1 

WRITE (6,1) K,HIST(T) , (LINE(M) ,M=1,J) 
110 CONTINUE 

1 FORMAT (IH ,I3,2X,F8.1, 1X,101A1) 
RETUR13 
END 


SUBROUTINE IlSTSTF (LUN ,HIST ,ERFLG , IFTC) 


THIS SUBROUTINE PUTS THE HlSTOGR/^1 ’ IIST’ INTO THE 
LABEL RECORD OF THE PICTURE IN 'LUN'. IT ALSO SETS 
THE FLAG THAT INDICATES THAT IT IS THERE. 'I PIC 

IS USED FOR AN I/O BUFFER. 'ERFLC IS SET TO ZERO 
IF THERE ARE NO ERRORS. 

SUBROUTINES USED - NONE 


REALM HIST (2 5 6) . FNUM 

INTEGER*2 ERFLG , IPIC (2048) , INUM(2) 

EQUIVALENCE (FKUM, INUM(l) ) 

ERFLG=0 

READ (LUN’1,EKD-110.ERR-120) (IPIC(I) . 1=1 , 102^) 

IPIC(512)=-1 

DO 100 1=2,512,2 

K=I/2 

FNUM=HIST(K) 

IPIC(512FI-1)=INUM(1) 

IPIC(512+I)=INUM(2) 

100 CONTINUE 

WRITE (LUN'1.END=110,ERR=120) (IPIC (I) , 1=1 , 1024) 
RETURN 

110 WRITE (7,1) LUN 
GO TO 130 

120 WRITE (7,2) LUN 

130 ERFLG- 1 .. 

1 FORMAT (’ END OF FILE ON UNIT *,13.’ IN ' ’HSTSTF* ' ) 

2 FORMAT (' I/O ERROR ON UNIT ',13,' IN "HSTSTF"') 
RETURN 

END 



onoonoooooooooooonnoo ooonrjnooooooo 


110 


SUBROUTINE INPUT (ST) 


THIS SUBROUTINE INPUTS THE DATA FROM THE FOUR 
POTENTIOMETERS. IT THEN MANIPULATES THEM SO THATt 
ST(1) IS A VALUE FOR M ON THE lUMGE (0,1); ST (2) 
IS THE LEFT SPREAD - IF IT IS POSITIVE, IT IS A 
DISPLACEMENT, NEGATIVE IS AN ANGLE (0,90); ST(3) 

IS THE RIGHT SPREAD - SAME EXPLANATION; ST(A) IS 
THE HEIGHT, H. 

SUBROUTINES USED - JOYSTK 


INTEGER*2 A,B,C,D 
REAL*4 ST (4) 

PAUSE 'TYPE C/R TO SET POTS’ 
CALL JOYSTK (A.B.C.D) 
ST(l)*A/4096.0-f0.5 
ST ( 2) = ( (9 0 . 0*B) / 2 048 . 0) 

IF (B.GT.O) ST(2)-B*0,01 
ST (3) -( (90. n-- 0 / 2048.0) 

IF (C.GT.O) ST(3)=C''0.0i 
ST(4)-1.0+(D/2048.0) 

IF (D.GT.O) ST (4) =1.0)0*0.01 

RETURN 

END 


SUBROUTINE INVCRV (TABLE, TEMP) 


THIS SUBROUTINE TAKES A TRANSFOPMATION 'TABLE', 
COMPUTES THE INVERSE OF THE ll'NS iOKMATION AND 
STORES IT RACK INTO 'TABLE', ULiHG 'TLMP' FOR 
WORK SPACE. 


THE INDEPENDENT VARIABLE IS ASSUME:) TO HAVE VALUES 
ON THE PJ^NGE (0,255). THE DEPENDENT VARIABLE IS 
ASSUMED TO HAVE THE SAME INTEGER IS\NGE. IF 
SEVERAL VALUES OF THE INDEPENDENT VARIABLE ARE 
MAPPED TO THE SAME VALUE OF DEPENDENT VARIABLE, 

THE INVERSE IS COMPUTED SUCH THAT THE DEPENDENT 


VALUE IS MAPPED BAXK TO ITS FIRST OCCURRENCE. IF 
THERE ARE NO OCCURRENCES OF A CERTAIN DEPENDENT 
VALUE IN THE INDEPENDENT VARIABLE, THE DEPENDI'.NT 
VALUE WILL MAP TO THE INDEPENDENT VARIABLE VniOSE 
DEPENDENT VALUE IS CLOSEST TO IT. 


SUBROUTINES USED - NONE 
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INTEGER*2 TABLE (256) , TEMP (256) 

>1 

K=TABLE(1)+1 
DO 100 I»1,K 
TEMP(I)-0 
100 CONTINUE 
M=K+l 

DO 140 I=M,256 
K=I“1 

110 IF (K. EQ. TABLE (J)) GO TO 130 
IF (K.LT.TADLE(J)) GO TO 120 
J-J+1 
GO TO 110 

120 IF ((K-TABLE(J-l)).LT.(TABLiJ(J)-K)) J”J-l 
130 TEMP(I)^J-1 
140 CONTINUE 

DO 150 1=1,256 
TABLE(I)-*TEMP(1) 

150 CONTINUE 
RETURN 
END 


SUBROUTINE JUG (IPIC J) 


THIS SUBROUTINE CONVERTS ’ITN’ ELEMENTS IN ROW 'J' 
OF ’IPIC' FROM BYTE DATA ON THE INTTIRVAL 
(-128, +127) TO INTEGER DATA ON . THE RANGE (0,255). 

SUBROUTINES USED - NONE 


INTEGER'’-2 IPIC (512, 4) 

DO 100 1=1, NW 

IF (IPIC(I,J) ,GE.128) IPIC(I,J)=IPIC(I,J)-256 
IPIC(I,J)=IPIC(I,J)+128 
100 CONTINUE 
RETURN 
END 


SUBROUTINE POTENT (HIST.ST.NREC) 


THIS SUBROUTINE CALLS ’INPUT' VmiCH GETS THE 
COEFFICIENTS 'ST' FROM THE POTENTIOMETERS, 

IF LS, RS AND H ARE ZERO, AND 1+0.5, IT GENER/\TES 
A UNIFORM HISTOGRAM AND RETUIGIS . IF LS , RS AND H 
ARE ZERO AND M IS NOT 0.5, IT GENEPJ^TES A 
HISTOGRAM THAT IS FLAT ON AN INTERVAL SYMMI^TRIC 
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ABOUT K, ZERO ELSEWHERE. IN ANY OTHER CASE, IT 
COMPUTES THE BREAKPOINTS OF THE PIECE- WISE LINEAR 
HISTOGRAM. 'H' IS THE Y COORDINATE AT ’M' . IF 
•LS' AND/OR 'RS' ARE LESS THAN ZERO, THEY 
REPRESENT ANGLES. THE X AND Y COORDINATES ARE 
THEN COMPUTED USING THE TRIGONOMETRIC FORMUIJV; 
GIVEN ANGLES B' and C AND INCLUDED SIDE A, 

1) A'=180-B'-C , 

2) B*A*SIN<B’)/SIN(A'). 

3) C-A*SIN(C')/STN(A') . 

IF 'LS' AND/OR 'RS' ARE GREATER THAN ZERO, THE 
BREAKPOINTS ARE (1,LS) AND/OR (256, RS). 'CRVGr.rt' 
IS CALLED TO GENERATE THE CURVE, GIVEN THE 
BREAKPOINTS, THE HISTOGRAM 'HIST' IS THEN SCALED 
TO THE CORRECT NUMBER OF POINTS USING ' NREC ' . 

SUBROUTINES USED - CRVGEN, INPUT 


REAL*4 HIST ( 256) , ST ( 4 ) , M , LS , RS , U 
100 CALL INPUT (ST) 

M=ST(1) 

LS»ST(2) 

RS=ST(3) 

H-ST(4) 

BB«H 

IF ((RS.NE. 0) .OR. (LS.NE, 0) .OR. (H.NE, 0)) GO TO 150 
IF (M.EQ.0.5) GO TO 230 
IF (M.LT.0.5) GO TO 110 
FL-2.0*M-1.0 
FH-1 . 0 
GO TO 120 
110 FL=0.0 
FH=2 . 0*H 

120 DH-(NREC-2)*8.0/ (FH-FL) 

IL-255.0*FL+1.5 
IH=255.0*FH+1.5 
DO 130 I=-l,256 
HIST(I)=0.0 
130 CONTINUE 

DO 140 I^IL.IH 
HIST(I)-DH 
140 CONTINUE 
GO TO 200 

150 ANGIKC-0. 0174532925 
PIBY2=1.5V0796327 
IF (LS.GT.O) GO TO 160 
THETA-PI BY2+LS‘'>ANGINC 
PHI-ATAM2(1 .0,M) 

A-M*COS (THETA) -'SIN (Pill ) / SIN (THRTA+PHI) 

AA“M*SIN (THETA) -SIN (PHI ) /SIN (THETA+PHI ) 

GO TO 170 
160 A- 0.0 
AA^LS 
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170 IF (RS.GT.O) GO TO 180 
ALPHA«PIBY2+RS*ANGINC 
BETA“ATAN2 ( 1 . 0 , ( 1 . 0 -M) ) 

C» ( 1 . 0-M) *C0S (ALPHA) *S1N ( BETA) /SIN ( ALPHA+BETA) 
C-l.O-C 

CC« o . 0-M)*SIN (ALPHA) *S IN (BETA) / SIN (ALPHAH-BETA) 
GO TO 190 
180 C«1 . 0 
CC-RS 

190 IA=255 0*A+1.5 
IB=255.0*M''1.5 
IC=255.0*C+1.5 

CALL CRVGEN (HIST, IA,IB,IC,AA,BB,CC) 

200 T0*0.0 

DO 210 1=1,256 
TO=TO+HIST(l) 

210 CONTINUE 

FH= (NREC-2)*2048. O/TO 
DO 220 1=1,256 
HIST(I)=HIST(I)>FH 
220 CONTINUE 
GO TO 250 

230 FH=(NREC-2)*8.0 
DO 240 1=1,256 
HIST(I)»FH 
240 CONTINUE 
250 RETURN 
END 


SUBROUTINE READO (LUN.NREC.ERFLC.LPIC) 


THIS SUBROUTINE READS \ BT.OCK OF DATA FROM RECORD 
’NREC’ OF FILE 'LUN'. THE DATA IS BYTE DATA ON THE 
RANGE (-128, +127) . THIS IS CONVERTED TO INTEGEIR 
DATA ON THE RANGE (0,255) AND STORED IN ’LPIC. 

IF 'ERFLG’=0, THERE WAS NO ERROR. 

SUBROUTINES USED - NONE 


INTEGER*2 ERFLG,J 
LOGICAL*! LPIC(4096) ,T(2) 

EQUIVALENCE (J,T(1)) 

ERFLG=0 

READ (LUK’NREC,END=110,ERR=120) (LPIC(I) , I»1 , 4096, 2 ) 

DO 100 1=1,4096,2 

T(2)=0 

T(1)=LPIC(I) 

IF (J.GE.128) J=J-256 
J=J+128 
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LPIC(I4l)rO 

LP1C(1)-T(1) 

100 CONTINUE 
RETURN 

110 WRITE (7.1) LUN 

GO TO no 

120 WRITE (7,2) LUN 
130 ERFLG-0 

1 FORMAT (’ END OF FILE ON UNIT ',13,' IN "READO"') 

2 FORHM: (' I/O ERROR ON UNIT ',13,' IN "READO"') 
R!:TURN 

END 


SUBROUTINE KECNliM (LP IC , NREC , P ICNA>I , LUN . ^.RFLG) 


THIS SlUiROin ILF. READS PART OF THE LABEL RECORD OF 
THE PICTURE ON FITE 'LUN' INTO TH?: BUFFER 'LPIC. 

IT THEN ASCEiTAINS THE SI2E OF T!X PICTURE - S, M, 
OR I., ATJD DETER’ UNIX TilE NUMBER OF RECORDS IN THE 
Fil-E 'MkEC. IT ALSO RLTUmiS THE N/EU: Or THE 
PICTURE 'PICXAM’. 'ERFLG' INDi CATES AH ERROR IF IT 
IS NOT ZERO. 


SUBROUTINES USED - NONE 


LOGICAI^l LPIC(4096) ,I'ICNAM(8) , S ,M,L 
lKTEGER-’'2 ERFLG , NREC 
DATA L.M.L/'S' , 'M' . 'L'/ 

NREC=--0 
ERFL(^“ 0 

READ (LUN’ i.END--110.F,KR=l20) (LPIC(I) , I»1 , 2048) 

IF (LPIC(A) .EQ.E) NRFC-10 
IF (LP1C(4) .FQ.H) NREC--34 
IF (LPIC(4) .EQ.L) NREC- 130 
IF (NREC.EQ.O) VJRITL (7,1) 

DO 100 

PICNAM(1)^'LPIC(I-4-8) 

100 CONTINUE 
RETURN 

110 WRITE (7,2) LUN 
GO TO 130 

120 \^TTE (7,3) LUN 
130 ERFLOl 

1 rORNlAT (’ INCORRfXT PICTURE SIZE’) 

2 FORMAT (* END OF FILE ON UNIT ',13,' IN "RECNUM"') 

3 FORMAT (’ I/O ERROR ON UNIT ',13,' IN "RECNUM"') 
RETURN 

END 
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SUBROUTINE RMSERR (REFHST , ACTHST , ERR) 


THIS SUBROUTINE COMPUTES THE ROOT MEAN SQUARE 
ERROR BETVIEEN THE REFERENCE HISTOGRAM ’REFHST’ AND 
THE ACTUAL HISTOGRAM ’ACTHST’. THE ERROR IS 
COMPUTED WITH THE HISTOGRAMS SCALED SO THAT THE 
SUMMATION OVER EACH HISTOGRAM IS ONE. 

SUBROUTINES USED - NONE 


REAL *4 REFHST (256) , ACTHST (256) , ERR 
T=0.0 

DO 100 1=1,256 
T=T+REFHST(I) 

100 CONTINUE 
T1=0.0 

DO 110 1=1,256 

DELT= (REFHST (I) -ACTHST (I )) /T 
T1=T1+DELT*DELT 
110 CONTINUE 

ERR=SQRT(Tl/256.0) 

RETURN 

END 


SUBROUTINE RWLBLO (LUN1,LUN2,SYM,SYSDAT, LABEL, ERFLG) 


THIS SUBROUTINE READS THE LABEL RECORD OFF OF THE 
INPUT FILE ’LUNl’, SUBSTITUTES THE APPROPRIATE 
PROGRAM IDENTIFIER ‘SYM’ INTO THE PICTURE NAME, 

AND WRITES THE LABEL RECORD INTO THE OUTPUT FILE 
’LUN2’. IT ALSO CLEARS THE FLAG THAT DENOTES THAT 
THE HISTOGRAM OF THE PICTURE IS IN THE LABEL 
RECORD. IT UPDATES THE SYSTEM DATE WORD, AS 
SUPPLIED BY ’SYSDAT’. 'LABEL' IS USED AS AN I/O 
BUFFER. 'ERFLG' INDICATES THERE HAS BEEN AN ERROR 
IF IT IS NOT ZERO. 

SUBROUTINES USED - NONE 


LOGICAL*! LABEL(4096) , SYM, SYSDAT (2) 

READ (LUN1'1,END=100,ERR=110) (LABEL (I) ,1=1, 2048) 
LABEL (16) =SYM 
LABEL ( 19) =SYSDAT(1) 

LABEL(20)=SYSDAT(2) 

LABEL (1024)=0 

WRITE (LUN2'I,END=130,ERR=120) (LABEL ( I) , 1=1 , 2048) 
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READ (LUNl* 2, END-100, ERR-1 10) (LABEL(I) ,1-1, 2048) 
WRITE (LUN2’ 2, END-130, ERR-120) (LABEL(I) ,1-1.2048) 
RETURN 

100 WRITE (7.1) LUNl 
GO TO 140 

110 WRITE (7,2) LUNl 
60 TO 140 

120 WRITE (7,2) LUN2 
GO TO 140 

130 WRITE (7,1) LUN2 
140 ERFLG-1 

1 FORMAT (’ END OF FILE ON UNIT ' ,13, ‘ IN ’ 'RWLBLO’ ’ 

2 FORMAT (’ I/O ERROR ON UNIT ’,13,’ IN "RWLBLO”’) 
RETURN 

END 


) 


SUBROUTINE STAT (ST, HIST) 


THIS SUBiiUTINE COMPUTES THE STATISTICS OF THE 
PICTURE ^OSE HISTOGRAM IS STORED IN ’HIST’ . THE 
STATISTICS COMPUTED INCLUDE THE MEAN, MODE, MEDIAN, 
MINIMUM AND MAXIMUM VALUES PRESENT, AND 
COEFFICIENTS OF SKEWNESS AND KURTOSIS. THE VALUE 
OF THE HISTOGRAM AT THE MODE IS ALSO COMPUTED. 
THESE VALUES ARE STORED IN ARRAY 'ST' AS WELL AS 
BEING OUTPUT TO THE TERMINAL DEVICE. 

SUBROUTINES USED - NONE 


REAL*4 HIST(256) ,ST(11) 

REAL*8 T1,T2,T3,T4 

TO-0 . 0 

Tl-O.O 

T2-0.0 

T3-0.0 

T4-0.0 

TM-O.O 

ST(7)-'’55.0 

ST(8)-0.0 

DO 100 1-1,256 

IF (HIST(I) .EQ.0.0) GO TO 100 
T— I**l 

T0»T0+HIST(I) 

T1»T1+HIST(I)*T 
T2-T2+HIST(I)*T*T 
T3-T3+HIST ( I ) *T*T*T 
T4-T4+HIST (I) *T*T*T*T 
TM=AMAX1 (TM , HIST (I) ) 
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IF (TM.EQ.HIST(D) ST(6)-T 
ST(8)-AMIN1(ST(8) ,T) 

ST(9)-AMAX1(ST(9) ,T) 

100 CONTINUE 
ST(5)-TM 
T1«T1/T0 
T2-T2/T0 
T3-T3/T0 
T4-T4/T0 
ST(1)~T1 
ST(2)»T2-T1*T1 

STnS-T3-3 0*Tl*T2+2 0*T1*T1*T1 

ST(4)»T4-4! 0*Tl*T3+6 ! 0*Tl*Tl*T2-3 . 0*T1*T1*T1*T1 

ST(10)«ST(3)/(ST(2)**1.5) 

ST(U)-(ST(4)/(ST(2)*ST(2)))-3.0 

Tl-0.0 

TM-TO/2.0 

11- ST(8)+1 

12- ST(9)+1 

DO 110 I«I1,I2 
T»»I-1 

Tl-TI+HIST(i) 

IF (Tl.GE.TM) GO TO 120 
110 CONTINUE 
120 ST(7)=T 

WRITE (7,1) ST 

1 FORMAT ('Ml - ‘.Ell. 4,' M2 = ' .Ell. 4. ’ M3 = '.Ell. 4, 
1’ M4 - ',£11.4./.' MODE » '.F8.1,' AT ‘.FS.l, 

2' MEDIAN » ’,F5.1,’ MIN » '.F5.1.‘ MAX « ’.F5.1./. 

3* A3 » '.E11.4.’ A4 - '.Ell. 4) 

RETURN 

END 


SUBROUTINE TRHIST .<TAB,HIST1,HIST2) 


THIS SUBROUTINE TRANSFORMS AN INPUT HISTOGRAM 
'HISTl' INTO AN OUTPUT HISTOGRAM 'HIST2* USING THE 
TRANSimMATION 'TAB'. 

SUBROUTINES USED - NONE 


INTEGER*2 TAB (256) 

REAL*4 HISTl (256) .HIST2 (256) 
DO 100 I»l,256 
HIST2(I)»0.0 
100 CONTINUE 

DO no 1=1,256 
J*TAB(I)+1 


oooooooooooo 00000000000 
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HIST2 ( J) -HIST2 ( J)+HIST1 (I) 
110 CONTINUE 
RETURN 
END 


SUBROUTINE UNIRMS (NREC, HIST. ERR) 


THIS SUBROUTINE COMPUTES THE ROOT MEAN SQUARE ERROR 
BETWEEN 'HIST' AND A UNIFORM HISTOGRAM WITH THE 
SAME NUMBER OF PIXELS AS DETERMINED BY ’ NREC ' . THE 
ERROR IS COMPUTED WITH THE HISTOGRAMS SCALED SO 
THAT THEY SUM TO 1. 'ERR' IS THE RESULTING ERROR. 

SUBROUTINES USED - NONE 


REAL*4 HIST (256) ,ERR 
PN’-(NREC-2)*8.0 
PNH= (NREC- 2) *2048 . 0 
T-0.0 

DO 100 I«l,256 
DELT= (PN-HIST (I) ) /PNH 
T-T-f'DELT*DELT 
100 CONTINUE 

ERR“SQRT(T/ 256.0) 

WRITE (7,1) ERR 

1 FORMAT ( * RMS ERROR WITH RESPECT TO UNIFORM - ’ , 

mi. 4) 

RETURN 

END 


SUBROUTINE WRITEO (LUN,NREC,LPIC,ERFLG) 


THE INPUT TO THIS SUBROUTINE IS AN ARRAY OF INTEGER 
DATA IN 'LPIC THAT HAS VALUES ON THE RANGE 
(0,255). THIS DATA IS THEN CONVERTED TO BYTE DATA 
ON THE RANGE (-128 ,+127) AND WRITTEN TO RECORD 
•NREC* IN FILE ’LUN*. 'ERFLG' INDICATES AN ERROR 
HAS OCCURRED IF IT IS NOT ZERO. 

SUBROUTINES USED - NONE . 


INTEGER+2 J, ERFLG 
LOGICAL*! LPIC(4096) ,T(2) 
EQUIVALENCE (J.T(l)) 


oooooooooooo 
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ERFLG-0 

T(2)-0 

DO 100 1-1,4096.2 
T(1)-LPIC(I) 

J— J-128 

IF (J.LT.O) J-J+256 
LPIC(I)»T(1) 

100 CONTINUE 

WRITE (LUN'NREC, END-110, ERR-120) (LPIC (I) ,1-1,4096,2) 
RETURN 

110 WRITE (7,1) LUN 
GO TO 130 

120 WRITE (7,2) LUN 
130 ERFLG-1 

1 FORMAT (’ END OF FILE ON UNIT M3,’ IN "WRITEO"') 

2 FORMAT (’ I/O ERROR ON UNIT ’,13,' IN ” WRITEO’") 
RETURN 

END 


SUBROUTINE XLATE (NREC , LUNl , LUN2 , ERFLG , IPIC , TABLE) 


THIS SUBROUTINE READS THE DATA FROM THE INPUT FILE 
’LUNl’, TRANSFORMS IT ACCORDING TO ’TABLE’, AND 
WRITES IT INTO THE OUTPUT FILE ’LUN2’. THIS IS 
DONE FOR THE ITHOLE PICTURE ACCORDING TO ’NREC’. 
’IPIC’ IS USED AS AN I/O BUFFER. ’ERFLG’ IS THE 
FLAG THAT SAYS THAT THERE WAS AN ERROR IF NOT ZERO. 

SUBROUTINES USED - READO, WRITEO 


INTEGER*2 IPIC (2048) .TABLE (256) .ERFLG 

DO 110 1=3, NREC 

CALL READO (LUNl , I , ERFLG , IPIC) 

IF (ERFLG. NE.O) GO TO 999 
DO 100 J-1,2048 
M-IPIC(J)+1 
IPIC (J) -TABLE (M) 

100 CONTINUE 

CALL WRITEO (LUN2, I, IPIC, ERFLG) 

IF (ERFLG. NE.O) GO TO 999 
110 CONTINUE 
999 RETURN 
END 


.TITLE CURSOR SAMPLER 

THIS SUBROUTINE SAMPLES THE JOYSTICK ON THE DIGITIZER 
AND RETURNS THE X, Y, AND Z COORDINATES OF THAT POINT. 
IT IS CALLED FROM FORTRAN IN THE FOLLOWING MANNER; 

CALL CURSOR (IX,IY,IZ) 

IX IS THE X COORDINATE OF THE JOY STICK CURSOR ‘ 
lY IS THE Y COORDINATE OF THE JOY STICK CURSOR 
IZ IS THE Z VALUE AT (IX, lY) AND IS OPTIONAL 


DEFINITIONS 


CREG 

-164120 

CONTROL AND STATUS 

XREG 

-164122 

SAMPLE X 

YREG 

-164124 

SAMPLE Y 

ZREG 

-164126 

Z OUTPUT 

JOXREG 

-164130 

JOY STICK X 

JOYREG 

-164132 

JOY STICK Y 

.GLOBL 

CURSOR 


.MCALL 
..V2. . 
.REGDEF 

.REGDEF, . .V2. . 



CURSOR; MOV 

JOXREG. @2 (R5) 

SAVE X 

MOV 

@#JOYREG,(a4(R5) 

SAVE Y 

MOV 

@2(R5) ,@#XREG 

SAMPLE X 

MOV 

@4(R5) ,(a#YREG 

SAMPLE Y 

CMP 

(R5).#2 

SAMPLE? 

BEQ 

OUT 

NO 

MOV 

#1,@#CREG 

GO 

LPl; BIT 

#200,(ay^CREG 

WAIT 

BEQ 

LPl 


MOV 

@#ZREG,@6(R5) 

GET Z 

OUT; RTS 

PC 


.END 



.TITLE 

JOYSTK SIMULATOR 



THIS SUBROUTINE SIMULATES A JOYSTICK WITH FOUR DEGRI 
OF FREEDOM, IT USES THE FOUR POTENTIOMETERS ON THE 


LABORATORY PERIPHERAL SYSTEM UNIT 
FOLLOWING MANNER: 

CALL JOYSTK (IA,IB,IC,ID) 


IT IS CALLED IN THE 


WHERE lA, IB, IC AND IB ARE INTEGER*2 VARIABLE THAT 
REPRESENT THE FOUR PARAMETERS. THE RANGE OF VALUES IS 
(-2048 ,+2047) . 
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WHEN THE SUBROUTINE IS CALLED, ALL OF THE LED'S ON THE 
LPS UNIT ARE BLANKED. THEN THE ROUTINE GOES INTO A LOOP 
THAT OPERATES IN THE FOLLOWING MANNER: 

1) BEGIN WITH CHANNEL 0 

2) SAMPLE CHANNEL, OFFSET IT TO (-2048.+20A7) 

3) MOVE BLANK OR MINUS SIGN INTO LEDy^2 AS REQUIRED 

4) COMPUTE INTEGER PERCENT (-99. +99) 

5) OUTPUT DIGITS TO LED'S 0 AND 1 

6) LOOK FOR INPUT FROM TTY -- IF NONE, GO TO 2 

7) SAMPLE CHANNEL, STORE IN CORRECT LOCATION 

8) INCREMENT CHANNEL # — IF NOT FINISHED, GO TO 2 

9) EXIT 


DEFINITIONS 


JSW 

-44 

; RT-11 JOB STATUS WORD 

ADS 

*170400 

; LPS A-D STATUS REGISTER 

ADB 

=170402 

; LPS A-D BUFFER REGISTER 

MINUS 

-1015 

; MINUS SIGN FOR LED#2 

PLUS 

=1017 . 

; BLANK IN LED #2 


.GLOBL 

JOYSTK 


.MCALL 
. .V2. . 

. .V2. . , .REGDEF, .TTINR, .TTYOUT 

• 

.REGDEF 





JOYSTK; 

CLR 

(3#ADS 


MOV 

#17, R1 


MOV 

#6,R2 

1$: 

MOV 

R1,@#ADB 

ADD 

#400. R1 


SOB 

R2,l$ 


TST 

(R5)+ 


BIS 

#10100, (ay^jsw 


MOV 

#4,R4 

2$; 

INC 

@#ADS 


TSTB 

@#ADS 


BPL 

. -4 


MOV 

(§#ADB.R2 


SUB 

#4000, R2 


MOV 

#PLUS ,@#ADB 


TST 

R2 


BPL 

3$ 


MOV 

#MINUS ,(a#ADB 


NEG 

R2 


CLEAR STATUS REGISTER 
CODE FOR BLANK TO LED 0 
SIX LED'S 
OUTPUT BLANK 
INCREMENT LED NUMBER 
LOOP 

BUMP ARGUMENT POINTER 
SET TTY BITS IN JSW 
FOUR CHANNELS TO SAMPLE 
START CONVERSION 
WAIT FOR READY 

MOVE VALUE TO R2 
OFFSET TO (-2048, +2047) 
MOVE BLANK TO LED #2 
IS R2 NEGATIVE? 

NO 

YES — PUT OUT MINUS SIGN 
TAKE ABSOLUTE VALUE 


THIS SECTION OF CODE CONVERTS THE CONTENTS OF R2 WHICH IS 
ON THE INTERVAL (-2048 ,+2047) TO A PERCENTAGE (-99, +99). 
IT ROUNDS OFF THE RESULT RATHER THAN TRUNCATE IT . 

A PERCENTAGE IS USUALLY COMPUTED BY R2*100/2048. THIS 
SUBROUTINE INSTEAD MULTIPLIES BY 3200 AND THEN TAKES THE 
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HIGH ORDER 16 BITS OF THE ANSWER. THIS ALLEVIATES THE 
NECESSITY FOR A DIVIDE. 


MUL 

#3200. ,R2 

ROL 

R3 

ADC 

R2 

CMP 

R2,#100. 

BLT 

4$ 

MOV 

#99. ,R2 

MOV 

R2,R1 

CLR 

RO 

DIV 

#10. ,R0 

ADD 

#400, RO 

MOV 

R0,@#ADB 

MOV 

R1,@#ADB 

.TTINR 

BCS 

2$ 

.TTYOUT 

INC 

(ay^ADS 

TSTB 

@#ADS 

BPL 

.-4 

MOV 

@#ADB,@(R5) 

SUB 

#4000.@(R5)+ 

INCB 

(ay^DS+l 

SOB 

R4.2$ 

TTYOUT 

#15 

BIC 

#10100, @#JSW 

RTS 

PC 

.END 


MULTIPLY BY 3200 
GET MSB OF LOW 16 BITS 
ROUND OFF ANSWER 
COMPARE TO 100 
OK 

SET TO 99 IF GREATER 
MOV TO LOW 16 BITS OF 32 
CLEAR HIGH 16 BITS 
CONVERT TO DECIMAL DIGITS 
PUT HI ORDER IN LED #1 
OUTPUT LED #1 
OUTPUT LED #0 
LOOK FOR TTY INPUT 
BRANCH IF NONE 
OUTPUT WHAT WAS INPUT 
SAMPLE SAME CHANNEL 
WAIT ON READY 

STORE VALUE 
OFFSET VALUE 

INCREMENT CHANNEL NUMBER 
LOOP 

OUTPUT CARRIAGE RETURN 

RESET JSW 

RETURN 


